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ABSTRACT 


In  the  early  chapters  (1-3)  an  expository  review  is  given  of 
Classical  and  Relativistic  Kinetic  Theory  of  simple  gases,  where  the 
Relativistic  Kinetic  Theory,  as  mainly  developed  by  J.L.  Synge  ["The 
Relativistic  Gas",  North  Holland  Publishers]  and  W.  Israel  ["Relativistic 
Kinetic  Theory  of  a  Simple  Gas",  Journal  of  Math.  Phys.  Vol.  4,  No.  9], 
is  that  of  a  simple  self-gravitating  gas  of  neutral  particles  with  no 
other  external  forces  undergoing  elastic  binary  collisions  and  where  the 
"mean-free  path"  is  space-time.  In  the  early  part  of  Chapter  4  small 
deviations  from  equilibrium  are  discussed  for  the  simple  gas  along  with 
certain  phenomenological  considerations.  In  tha  latter  part  of  Chapter 
4  new  explicit  expressions  for  the  transport  coefficients  (i.e.  Conductiv 
ity,  Bulk  viscosity,  and  Shear  viscosity)  of  the  simple  relativistic  gas 
are  given.  In  Chapter  5  certain  consequences  of  the  existence  of  Bulk 
viscosity  and  its  asymptotic  temperature  dependence  are  briefly  discussed 
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CHAPTER  I 
Equilibrium  Stat e 


In  this  chapter  a  brief  discussion  of  the  statistical  nature 
of  the  equilibrium  configuration  of  a  simple  gas  is  discussed  both  for 
the  Classical  and  Relativistic  cases. 

Given  a  gas  consisting  of  n  particles  where  n  is  quite 
large,  we  would  like  to  describe  its  macroscopic  behavior  in  terms  of 
its  microscopic  variables.  To  present  a  precise  account  of  its  behavior 
would  clearly  be  a  gargantuan  task  because  of  the  enormity  of  the  initial 
data  and  variables  involved.  This  leads  one  to  a  statistical  approach  to 
the  problem.  In  this  section  the  usual  statistical  approach  used  in  the 
derivation  of  the  Maxwell-Bolt zmann  distribution  function  will  be 
discussed  for  the  Classical  and  Relativistic  cases. 

The  statistical  approach  takes  the  form  of  specifying  some 
distribution  function.  Instead  of  knowing  each  particle's  exact  behavior 
i.e.  its  position  and  momentum,  we  will  now  know,  via  the  distribution 
function,  that  a  given  number  of  particles  have,  for  example,  a  momentum 
within  a  certain  range  centered  about  some  value  of  momentum.  In  the 
classical  case  the  space  in  which  we  operate  in  a  6-dimensional  one.  It 
is  called  the  p-space.  Even  though  each  particle  has  6  degrees  of 
freedom  (excluding  spin  effects,  etc.),  and  there  are  n  of  these  parti¬ 
cles,  we  superimpose  them  all  one  one  6-dimensional  grid.  Mathematically 


2 


then  we  have  some  distribution  function  N(r,v,t)  which  is  a  function 

of  position,  velocity,  and  time.  Given  an  elemental  region  in  our 

3  3 

y-space,  i.e.  d  rd  v  (1)  about  some  point  (r,v)  ,  then 

3  3 

N(r,v,t)d  rd  v  represents  the  number  of  particles  that  have  coordi- 

3  3 

nate  values  that  fall  within  the  elemental  volume  d  rd  v  at  time  t  . 

It  is  in  this  setting  that  the  Boltzmann  Transport  equation  is  obtained. 

But  in  order  to  get  a  more  fundamental  insight  into  what  the  equilibrium 
configuration  means  statistically,  we  will  enlarge  our  space  from  a  6- 
dimensional  one  to  a  6n-dimensional  space.  This  new  space  is  called 
the  T-space.  In  the  y-space  representation  n  points  determine  the 
state  of  the  gas.  But  in  T  space  one  point  specifies  the  whole  system. 
Suppose  we  interchange  two  particles  from  two  cells.  This  would  still 
give  us  the  same  configuration  in  y-space,  i.e.,  the  same  N  ,  but  when 
this  interchange  is  viewed  in  T-space  it  leads  to  a  new  distinct  point. 
Because  of  this  we  see  that,  associated  with  a  given  distribution  function 
N  ,  or  configuration  in  y  space,  we  have  a  multiplicity,  or  a  cloud,  of 
points  in  T  space.  In  fact  this  multiplicity  should  be  expected,  because 
we  should  really  be  starting  in  the  T-space,  and  since  y-space  is  a  kind 
of  projection  of  T  space  down  onto  one  6  -dimensional  grid,  one  expects 
to  lose  information  in  the  process. 

The  equilibrium  state  for  a  gas  will  now  be  associated  with  that 
state  that  maximizes  the  "volume"  that  the  system  takes  up  in  T  space. 
What  we  are  saying  is  that,  given  a  system  with  certain  macroscopic  limita¬ 
tions  (e.g.  fixed  energy,  number,  etc.),  the  equilibrium  state  will  be  that 
which  allows  for  the  greatest  number  of  permutations  of  the  microscopic 
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variables  and  yet  leaves  the  macroscopic  system  unchanged  (2).  This  is 
equivalent  to  saying  that  our  system  will  have  the  most  probable  distri 
bution. 


The  problem  therefore  becomes  one  of  maximizing  a  volume  or 
maximizing  the  probability  associated  with  any  given  distribution.  The 
approach  that  will  be  taken  is  that  of  maximizing  the  probability  P  , 
where  P  stands  for  the  probability  of  having  a  given  distribution  . 
Thus  we  essentially  can  stay  in  y- space. 

Divide  up  y-space  into  elemental  cells  all  having  the  same 
magnitude  of  extension,  and  index  them  by  1 ,2  ,  • • •  , i  ,  •  • •  .  Let  us  now 
take  the  n  particles  and  distribute  them  out  amongst  the  cells,  where 
we  write  n_^  as  the  number  of  particles  in  the  ith  cell.  But  the 
particles  can’t  distribute  themselves  in  any  fashion,  they  must  first 
satisfy , 

(1-1)  l  n  =  n  . 

i 

Since  we  are  assuming  that  the  total  energy  is  fixed  then  we  must  also 
satisfy , 

(1-2)  l  n.£l  =  E  , 

i 

where  e.  is  the  energy  per  particle  in  the  ith  cell.  It  is  a  well 

i 

known  fact  from  combinatorics  that  the  probability  for  a  given  distribu 
tion  (3),  for  fixed  population  numbers  n^  ,  • • • ,n^  ,  • • • ,  is 


(1-3) 


T3  — 
r 


x  const. 


-  4  - 


_ nj _ 

n, !n0 ! • • *n. ! 

12  i 

Taking  the  log  of  equation  1-3,  we  get, 

(1-4)  log  P  =  log  (n!)  -  £  log  (n.!)  +  const. 

Assuming  that  n_^  are  all  quite  large  we  can  make  use  of  the  approxima¬ 
tion 


(1-5) 


log  (n!)  =  n  log  n  -  n  , 


which  follows  from  Sterling's  formula.  Substituting  equation  1-5  into 
equation  1-4,  and  noting  equation  1-1  we  get, 

(1-6)  log  P  =  n  log  n  -  £  n.  log  n.  +  const. 

i 

Since  the  most  probable  distribution  means  it  is  more  probable  relative 
to  any  other  distribution,  let  us  then  compare  any  two  distributions  of 
n  particles  to  obtain  their  relative  probability.  .  By  means  of  equation 
1-6  we  see  that, 


(1-7) 


log(P/P’>  =  -  l  n  log  n  +  l  n'  log  n' 

j  3  3  j  j  3 


where  n.  and  n!  are  the  occupation  numbers  of  the  two  distributions, 

J  J 

Letting 


Q  =  "  I  n  lo 
j 


g  n, 


Q'  = 


l  l°g 

j 


(1-8) 
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it  can  be  said  that  the  distribution  (n^, • • • ,n^ • • • }  is  more  probable 

than  the  distribution  (n! , • • *n!  ,  •  •  •  }  i.e.  P  >  P’  if  Q  >  O'  .  Simi- 

1  J 

larly,  P’>P  if  Q*  >  Q  ,  and  P=P*  if  Q'=Q.  In  comparing  two 
distributions  there  are  two  conditions  which  are  required.  Firstly,  the 
total  number  of  particles  must  be  the  same.  Secondly,  the  division  of 
y-space  into  elemental  cells  must  be  done  in  the  same  fashion  for  both 
distributions.  Since  we  have  divided  y-space  into  cells  of  equal  exten¬ 
sion  we  label  the  magnitude  of  extension  as  eq  ,  where  e  is  the  spatial 
volume,  q  the  velocity  (or  equivalently  the  momentum)  volume.  We  can 
thus  write 


(1-9)  n.  =  eq  N.  , 

J  3 

where  N.  is  the  value  of  the  distribution  function  over  the  i  th  cell. 
1 

Using  equation  1-9  the  expression  for  Q  can  be  rewritten  as 


(1-10) 


Q 


eq  )  N.  log  N. 
L  J  3 


eq  n  log  eq 


Going  back  to  equations  1-7,  1-8,  and  in  view  of  equation  1-10, 
it  can  be  seen  that  when  probabilities  are  compared,  we  are  essentially 
comparing  expressions  of  the  form 

(1-11)  F  =  -  l  N  (log  N  )  eq  . 

F  is  therefore  a  measure  of  the  relative  probability.  Demanding  the 


usual  conditions,  equation  1-11  can  be  rewritten  as 


' 


(1~12) 


F 


r 


-  6  - 


N  log  N  d^r  d^v  , 

3  3 

where  e  ,  and  r|  have  become  d  r  ,  and  d  v  respectively.  This 

quantity  F  ,  which  is  a  measure  of  the  chaos  of  the  system,  is  called 

3 

the  entropy  integral.  Considering  some  fixed  elemental  volume  d  r  , 
over  which  F  will  be  maximized,  the  equilibrium  condition  is 

(1-13)  5F  =  0  , 


where  the  previous  constraints  (equations  1-1, 1-2)  for  the  continuous 
case  are 


n  = 


3  3 

N  d  r  d  v 


(1-14) 


where 


6n  =  0 


|  3  3 

E  =  Nedrdv  where 


SE  =  0 


Varying  N  with  5N  arbitrary,  and  using  the  above  constraints, 
equation  1-14  becomes 


(1-15) 


(log  N  +  1)  6N  =  A6N  -I-  3e6N 


where  A  and  3  are  constants.  Since  the  6N  are  arbitrary  equation 
1-15  can  be  put  into  the  form 


N  =  C  e 


(1-16) 


- 
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m  3/2 

It  can  be  shown  that  C  =  n(— -)  and 

r\  £v  I 

p 

Recalling  that  e  =  —  >  we  see  that  equation  1-16 


Boltzmann  distribution. 


is  the 


(4). 

usual  Maxwell- 


In  the  Relativistic  case  the  equilibrium  condition  is  arrived 
at  by  essentially  the  same  methods  as  were  used  for  the  Classical  case. 
That  is,  the  equilibrium  condition  is  characterized  by  6F  =  0  ,  where 
F  has  the  same  form  as  the  classical  case.  But  in  the  Relativistic  case, 
the  formulation  of  y-space  needs  closer  scrutiny.  Also  this  approach 
must  make  manifest  its  covariance.  Before  the  relativistic  derivation  can 
be  given  there  are  certain  preliminaries  which  must  first  be  given.  We 
say  that  space-time  forms  a  4-dimensional  manifold.  Each  point  in  this 
manifold  is  called  an  event,  labelled  ,  where  the  global  space-time 

structure  need  not  be  Minkowskian.  Consider  a  particle,  and  at  an  event 
x^  associated  with  the  particle  is  the  4-momentum  P^  .  The  4-momentum 
of  the  particle  lies  in  the  "plane”  tangent  at  x^  to  space-time.  We 
call  this  plane  the  momentum  space.  The  momentum  space  is  Minkowskian. 

At  the  event  x'^  we  can  form  a  tetrad  in  Momentum  space.  The  vectors  in 
the  tetrad  are  tangent  to  the  coordinates  net  in  space-time  at  event  x^  . 
Since  all  of  our  work  will  be  carried  on  in  a  local  Minkowskian  frame  in 
space-time  we  see  then  that  each  vector  in  the  momentum  tetrad  is  parallel 
to  a  corresponding  vector  in  the  tetrad  we  pick  in  space-time. 

Of  all  those  vectors  that  exist  in  momentum  space,  the  only  ones 

of  interest  will  be  the  future  pointing,  time-like  vectors,  i.e. 

2  X 

P  P^1  =  -in  ,  and  P  f  >  0  .  Since  in  our  case  we  will  be  concerned  with 

y 


. 
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the  situation  where  all  the  particles  have  the  same  constant  proper 
mass  it  is  seen  that  P^Py  =  -m2  describes  a  pseudosphere  (5)  of 
radius  m  .  Thus  all  the  particles  in  the  system  have  a  4-momentum 
vector  which  terminates  on  the  pseudosphere  of  radius  m  .  Polar  co¬ 
ordinates  can  be  introduced: 


(1-17) 


P  =  m  sinh  x  sin  0  sin  <J> 
2 

P  =  m  sinh  x  sin  0  sin  4) 
3 

P  =  m  sinh  x  cos  0 
P^  =  m  cosh  x 


where  0_<x<o°  ,  O<_0<_tt  ,  0  <_  <j)  <  2tt  ,  (if  we  were  admitting 

variable  proper  masses  then  0  <  m  <  °°)  . 

In  the  classical  case  the  distribution  function  N  is  a 
function  of  r,v,  and  t  ,  and  N  is  defined  such  that 

v  =  N(r,v,t)  d2r  d^v 

represents  the  number  of  particles  an  observer  would  see  if  he  constructed 
3 

a  volume  d  r  about  a  point  r  and  counted  the  number  of  particles  whose 

3 

velocity  vectors  were  within  a  range  d  v  of  v  .  Because  of  the  fact 
that  in  Newtonian  Theory  "space”  degenerates  into  a  plane,  our  volumes 

have  an  absolute  character.  But  this  doesn’t  quite  work  out  in  Relativity. 

Depending  on  how  the  time-like  unit  vector  is  chosen,  there  will  be  many 
3  flats  orthogonal  to  the  normal.  It  should  be  pointed  out  that  momentum 
space  will  be  used  instead  of  velocity  space  (6). 


p 
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The  volume  that  replaces  the  classical  volume  d  v  ,  will  be 
a  3-cell  lying  on  the  pseudosphere,  so  that  given  the  4-momentum  P^ 
of  a  particle  it  will  terminate  on  the  pseudosphere  and  all  the  parti¬ 
cles  with  4-momemta  that  lie  within  a  given  range  of  P^  will  terminate 
on  an  element  of  surface  of  the  pseudosphere.  Pictorically  we  have 
Figure  1,  where  A  is  a  3-cell  about  a  point  on  the  pseudosphere. 


To  see  how  the  distribution  function  is  defined  within  the 

Relativistic  framework,  let  us  consider  at  some  event  a  time-like 

unit  vector  n  and  orthogonal  to  n  an  elemental  3-flat  ds  which 
p  y 

will  serve  as  our  target.  We  want  to  know  how  many  particles  there  are 
with  the  properties  that  their  4-momenta  lie  within  a  given  elemental 
range  centered  about  some  given  momentum  P^  and  that  they  cross  ds  . 
To  take  into  account  the  fact  that  n^  is  not  coincident  with  P1-*  ,  and 
hence  our  3-flat  ds  is  not  normal  to  P^  ,  we  project  the  elemental  3- 
area  on  the  pseudosphere  along  n^  down  onto  a  3-flat  parallel  to  ds  . 
Let  this  projected  3-volume  be  called  dQ  .  Pictorically  we  have  figure 
2. 


sphere , 


Noting  that  (— )  (— )  =  -1 

mm 

then  by  differentiating  we  get 


describes 


the  surface  of  the  pseudo- 


)  d( 


=  0  . 


P^1 

Therefore  ( — )  is  orthogonal  to  the  tangent  on  the  pseudosphere, 
m 
ph 

hence  ( — )  is  orthogonal  to  the  pseudosphere.  Denoting  the  cell  A  by 


mdm  ,  then  dW  and  doo  can  be  related  by  means  of  the  Projection  theorem 


. 

■ 
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pseudosphere 


FIGURE  1 


FIGURE  2 
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(7) .  That  is , 


(1-18) 

dft  = 

=  mdco 

y  1 

'  m  ny  * 

or 

d^  = 

=  da) 

Pyn  | 

y 

To  relate 

dfi  and 

dco 

in  terms  of 

our  polar 

coordinates 

(i .  e . 

equation  1-17) 

,  we  use 

n 

y 

as  our  axis 

off  itfhich 

X  is  measured 

(i  .e. 

X  =  0  along 

n  ) . 

y 

Since  n 

y 

lies  al 

ong  the  axis ,  thus 

nl  =  n2 

=  n^  =  0  and 

n4  =  " 

1  .  Equation 

1-18  therefore 

becomes 

dfi  = 

du)  |  - 

m  cosh  x i 

or  , 

(1-19) 


dfi  =  m  cosh  x  dm 


In  the  3-flat  orthogonal  to  n  we  are  just  using  normal 

H 


spherical  coordinates,  where  the  radius  is 


(1-20) 


R  =  m  sinh  X 


An  element  of  volume  in  spherical  coordinates  is 


(1-21) 


dS!  =  R  sin  9  d  R  d0  d<}> 


Substituting  equation  1-20  into  equation  1-21  we  get 


3  2 

dSl  =  m  sinh  x  cosh  x  sin  9  dx  d9  d(}> 


(1-22) 
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Therefore  by  equation  1-19  we  get, 

(1-23)  dm  =  m^  sinh^x  sin  6  d  x  d6  dcj) 

The  quantity  dm  is  called  the  absolute  2-content  of  a  3-cell 

on  the  pseudosphere.  The  terminology  "2-content"  follows  from  the  fact 

2 

that  dm  has  dimensions  of  (Mass)  ,  which  should  not  seem  strange 

since  we  identified  the  volume  of  a  3-cell  as  mdco  which  we  know  has 

3 

dimensions  (Mass) 

For  the  sake  of  future  reference  the  relationship  of  dm  to 
the  elementary  solid  angle  in  momentum  space  will  be  discussed.  This 
will  be  of  use  when  discussing  the  problem  of  scattering,  where  we  want 
to  know  how  many  particles  are  scattered  into  a  given  solid  angle  (i.e. 
a  given  momentum  range) .  Normally  a  solid  angle  is  defined  to  be  equal 
in  magnitude  to  the  surface  area  it  subtends  on  a  unit  sphere.  Similarly, 
given  a  unit  pseudosphere,  i.e.  m  =  1  ,  then  by  equation  1-23  we  get 

2 

(1-24)  dm  =  sinh  x  sin  9  dx  d9  dcj> 

2 

So  sinh  x  sin  9  dx  d9  d<£  becomes  our  elementary  solid  angle. 

Also  dm  now  represents  the  numerical  value  of  the  volume  of  a  3-cell  on 

the  pseudosphere  since  mdu)  =  dm  .  If  we  don't  have  a  unit  mass  then 

the  elementary  solid  angle  is  .  Many  authors  (8)  make  use  of  the 

m 

simplicity  of  equation  1-24  by  dealing  with  the  normalized  momentum  vectors, 
u 

PH  u 

i.e.  ( — )  ,  instead  of  P  .  This  has  the  effect  of  operating  on  the  unit 

m 


' 

■- 

■ 

■ 
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pseudosphere . 

Getting  back  to  the  distribution  function  N  ,  we  define  it 

such  that 


v  =  N(x,P)  ds  d^  , 


represents  the  number  of  world  lines  that  intersect  ds  with  normal 

n  and  such  that  the  4-momenta  fall  within  the  3-cell  dQ  ,  where 
h  5 

dfi  =  In  |  dco  .  da)  is  an  invariant  but  dG  dpends  on  the  choice  of 
n ^  .  It  turns  out  that  N(x,P)  is  an  invariant  i.e.  it  does  not  depend 
on  which  3-flat  we  take  our  sample  populations  from.  The  proof  of  the 
invariance  of  N  rests  essentially  on  two  points.  The  first,  which  has 
already  been  mentioned,  is  the  invariance  of  dm  .  The  second  point  is 


that  if  we  look  at  the  same  configuration  of  particles  with  respect  to 

What  we  are 

saying  is  that  the  product  ds  dfi  is  also  an  invariant.  The  reason  for 


two  different  normals  n  and  n*  ,  then  ds  dfi  =  ds’d^' 

y  y 


stipulating  a  fixed  configuration  is  that  if  our  comparison  of  two  distri¬ 
bution  functions  is  to  be  meaningful,  then  we  should  be  referring  to  the 
same  set  of  particles.  This  demand  will  be  guaranteed  by  the  way  we  con¬ 
struct  our  elemental  3-flats  ds  d^  ,  ds '  and  d£2'  . 


Consider  a  local  tube  of  world  lines  characterized  by  momentum 
Pli  at  event  x^  ,  (see  figure  3) .  Let  n^  be  some  time-like  unit 
vector  at  event  x^  ,  and  let  ds  be  a  cross  section  of  the  world-tube 
normal  to  n  ,  (see  figure  4) ,  and  suppose  v  is  such  that 
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FIGURE  3 


FIGURE  4 
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(1-25) 


v  =  N  (x  ,P  )  ds  dH 


Consider  another  time-like  unit  vector  n’  at  event  x^  , 

y 

and  let  ds ’  be  a  cross  section  of  this  same  tube  normal  to  n'  then 

y 

suppose  we  have 


(1-26) 


v'  =  N' (x,P)  ds'  dQ' 


We  can  now  guarantee  that  we  have  the  same  configuration  of 
particles  if  we  say  that  dfi '  and  dfi  correspond  to  the  same  2-content 
dco  .  This  is  we  can  start  with  dft  and  obtain  our  dfif  by  projecting 


dfi  up  to  doo  along  n  and  then  projecting  down  along  n 

r* 


to  obtain 

y 


d$7'  .  Pictorically ,  we  have  figure  5.  By  the  projection  theorem  we  have 


(1-27) 


ds '  I P^n 'I  =  ds  j  P^n  I  } 

y  1  y1 


and  we  also  have  by  equation  1-18 

dtt' 


(1-28) 


dco  = 


dQ, 


|pV| 

|pyn  | 

1  y ' 

1  y 1 

Equations  1-27,  1-28  yield  the  result 


(1-29)  ds  did!  =  ds '  dO*  . 

Since  we  are  dealing  with  the  same  particle,  v  =  v*  ,  so 
that  using  equation  1-29,  equations  1-25,  1-26  yield  N  -  N ’  .  This 
tells  us  that  the  distribution  function  N  is  independent  of  our  choice 


-  16  - 


m 


FIGURE  5 
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of  n 

U 


As  was  mentioned  earlier  the  derivation  of  the  equilibrium 
distribution  function  is  analogous  to  the  classical  case.  But  the 
cells  e  and  q  now  represent  elemental  cells  in  the  3-flats  of 
space-time  and  momentum  space  respectively,  where  they  are  normal  to 
some  given  timelike  unit  vector,  i.e. 


By  arguments  similar  to  those  used  for  the  classical  case, 
we  will  maximize  the  following  entropy  integral 


(1-30) 


F 


N  log  N  ds  dQ  , 


and  the  condition  for  equilibrium  will  be 

(1-31)  6F  =  0  . 

The  only  thing  left  is  to  state  the  constraints.  In  the  classi¬ 
cal  case  we  could  vary  N  as  long  as  in  so  doing  we  didn't  violate  the 
conditions  of  constant  number,  and  energy.  As  it  turns  out  the  relativ¬ 
istic  constraints  are  of  a  local  nature.  That  is,  the  fixed  quantities 
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(or  predeterminations)  that  must  be  conserved  when  varying  N  are 
specified  locally.  To  establish  what  these  predeterminations  might  be 
we  make  use  of  a  rather  ingenious  construction  due  to  J.L.  Synge  (9). 

The  key  to  this  construction  is  that  it  utilizes  the  time-like  character 
of  the  4-momenta  of  the  particles.  We  pick  our  target  by  considering 
at  some  event  x^  a  future  pointing  null  cone,  and  we  consider  a  cross 
section,  s  ,  of  the  null  cone  which  is  normal  to  some  time-like  unit 
vector  n^  .  Let  J  be  the  region  of  the  null  cone  between  the  event 
x^  and  the  cross  section  s  . 


The  beauty  of  this  construction  is  that  the  world  line  of  a 
particle  crossing  £  must  then  cross  s  ,  and,  as  long  as  we  have  con¬ 
servation  of  numbers  in  a  collision,  then  collisions  in  the  space-time 
volume  bounded  by  £  ,  and  s  will  not  alter  this  fact,  i.e.  world 
lines  cannot  leak  out  of  this  volume  (see  figure  6).  Because  of  the 
conservation  of  4-momenta  at  each  collision  we  see  that,  given  a  momen¬ 
tum  flux  across  £  ,  then  s  must  have  the  same  momentum  flux. 


Our  constraints  therefore  take  the  form  guaranteeing  the  rela¬ 
tionship  between  the  world  lines,  and  the  momenta  crossing  £  ,  and  s  . 
Everything  that  has  been  said  up  to  this  point  is  not  necessarily  of  a 
local  nature.  But  now  we  want  to  consider  an  infinitesimal  target  ds  so 
that  s  must  be  close  to  x^  .  In  fact  we  assume  that  the  distribution 
function  N  ,  the  numerical  flux  vector  N  dco  ,  and  the  energy- 


momentum  tensor  T 


pv 


P^PVN  do)  are  smooth  continuous  functions  so  that 


their  values  on  ds  are  approximately  those  at  x 


P 


It  is  at  this  stage 
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n 


FIGURE  6 
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that  the  argument  becomes  one  of  a  local  nature,  since  N  ,  Mm  ,  and 
T^V  are  evaluated  at  x^ 


Suppose  that  the  number  of  world  lines  crossing  £  is  s 


ome 


fixed  quantity  v  .  Then  this  must  also  be  the  same  number  crossing 

ds  .  Since  the  number  of  world  lines  crossing  ds  is  -M  n  ds  then 

U 

we  can  write 


(1-32) 


v  =  -  M^n  ds 

y 


Since  M^1  =  P1J  N  dco  ,  equation  1-32  becomes 


(1-33) 


v  =  ds 


-  PU  n  N  dm 
U 


But  P^  and  n^  are  both  time- like  which  means  that  -P^  n  >  0  ,  so 

y 

that  equation  1-33  can  be  rewritten  as 


v  =  ds  N  |P,J  n  |  du 


By  equation  1-18  the  above  becomes 


(1-34) 


v  =  ds 


N  dQ. 


Now  suppose  that  the  momentum  flux  crossing  £  is  some  fixed 


quantity,  say  M 


^  .  But  the  momentum  flux  across  ds  is  (~T^ V  n^  ds) 


Since  T 


nyv 


P^1  PV  N  dm  ,  we  can  write 


.  ' 
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(1-35) 


uv 

-  T‘  n  ds  =  ds 
v 


u  I  V  I 

P  P  n  dm  =  ds 
1  v 1 


P1J  N  dfi  . 


Therefore  we  can  write, 


(1-36) 


My  =  ds 


N  Py  dO, 


Thus  the  variations  of  N  must  be  such  that  they  are  constraints  (or 
predeterminations),  equations  1-34,  1-35  are  not  violated. 

Using  lagrangian  multipliers ,  the  condition  6F  =  0  reads 


(1-37) 


(log  N  +  1)6N  =  c  6 N  +  3y  Py6N 


where  c  and  3^  are  independent  of  P^  ,  but  are  in  general  functions 
of  .  Equation  1-39  becomes 


(log  N  -  c  -  3^  PP  +  1) 6N  =  0 


which  implies  that 


i  3  P 

IT  C-l  P 

N  =  e  e 


y 


C —  1 

Letting  e  =  a(x)  be  finally  get 


(1-38) 


N(x,p)  =  a(x)  e 


6  P 
y 


y 


where  8  must  be  time-like  to  guarantee  convergence.  Using  equation 

y 

1-38  it  can  be  shown  that  we  get  the  numerical  flux,  and  energy  momentum 
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tensor  of  a  perfect  fluid  with  the  hydrodyncmic  4-velocity  properly 
defined.  But  we  save  this  for  a  later  section  where  the  conservation 
laws  will  be  derived  in  a  natural  way  from  the  Boltzmann  equation. 


23 


CHAPTER  II 

Boltzmann  Collision  Equation 


In  this  chapter  the  usual  derivation  of  the  Boltzmann  Colli¬ 
sion  Equation  is  given  for  both  the  Classical  and  Relativistic  Cases. 

In  the  previous  chapter  we  neglected  the  more  general  case 
of  non-equilibrium.  If  we  have  a  non-equilibrium  situation  we  know  it 
will  eventually  reach  equilibrium,  (i.e.  the  distribution  function  will 
approach  the  Maxwell-Bolt zmann  distribution) ,  but  to  this  point  we  can¬ 
not  describe  what  the  behavior  of  the  system  will  be.  For  example  it  is 
in  this  approach  to  equilibrium  where  the  transport  effects  occur  as  a 
means  for  smoothing  out  the  non-uniformities  responsible  for  the  non¬ 
equilibrium  situation.  It  is  at  this  stage  that  the  Boltzmann  transport 
equation  is  introduced  to  explain  the  more  general  behavior  of  the  gas. 

The  Boltzmann  transport  equation  describes  the  evolution  of  the  distribu¬ 
tion  function  by  taking  into  account  the  effect  collisions  have  on  N  . 

To  make  this  clearer,  we  will  give  the  following  rather  crude  argument 
which  will  serve  both  to  illustrate  the  above  statement,  and  to  introduce 
the  Boltzmann  transport  equation. 

Suppose  that  at  r,v  we  have  some  hypothetical  observer  at 

3  3 

time  t  .  About  him  we  have  the  volume  d  r  d  v  .  Suppose  he  observes 
some  density  distribution  N  of  particles  in  this  volume.  Further  suppose 
\<re  stipulate  that  there  are  no  collisions.  We  can  view  the  time  evolution 


. 


' 

. 
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of  these  points  as  a  mapping  of  y-space  into  y-space.  After  an  infini- 
tesmal  time  dt  our  "observer"  has  moved  with  the  volume  from  r,v  at 

'  r-*j 

t ime  t  to 


r  + 


dr 

nsj 

dt 


dt 


v  + 


dv 

r*>j 

dt 


dt 


5 


3  3 

at  a  time  t  +  dt  ,  and  the  volume  has  become  d  r'  d  v*  .  But  by 

3  3  3  3 

Liouvilles  theorem  drdv=drTdv'  .  Since  we  are  talking  about 
the  same  particles,  we  see  that,  under  the  assumption  of  no  collisions, 
our  observer  will  see  no  change  in  the  density  of  points.  That  is, 


(2-1) 


dN 

dt 


0 


,  where 


is  the  total  derivative. 


Since  N  is  a  function  of  r,v  ,  and  t  ,  we  write  equation 

2-1  as, 


dN  _ 

3N 

3t  + 

3N 

dr 

dt  + 

3N 

dv 

dt 

3r  * 

r^j 

3  v 

d7 

(2-2) 

m 

3N 

37 

P  • 

3N  , 

37  + 

F  • 

3N  . 
3v 

n+j  /-*w* 


What  we  are  saying  is  that,  under  the  condition  of  no  collisions,  the 
density  per  unit  "volume"  of  particles  or  points  in  y-space,  behaves  like 
an  incompressible  fluid.  But  if  we  allow  for  close  binary  collisions 
then  there  will  be  particles  which  will  leak  out  our  volume  due  to  these 
collisions.  If  some  leak  out  of  one  elemental  volume  they  then  must  leak 
into  another.  Therefore,  in  general,  our  observer  moving  with  the  volume 
will  not  find  equation  2-2  satisfied,  but  instead  he  will  observe  a  change 


% 


• 
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in  N  due  to  a  net  change  in  the  number  of  particles  in  the  elemental 
volume.  To  take  this  into  account  we  modify  equation  2-2  by  adding  the 
customary  phenomenological  collision  term,  i.e.  for  collisions  equation 
2-2  becomes 

fry  ^  3N  9N  ^  3N  T 

(2-3)  — -+v*  —  +  —  •—  =J  .  -  J  , 

3t  ~  9r  m  3v  gam  loss 

where  J  .  -  J,  gives  us  the  net  gain  or  loss  of  particles  per  unit 

gain  loss 

volume  in  y-space.  Equation  2-3  is  called  the  Boltzmann  Transport  equa¬ 
tion.  The  phenomenological  tern  J  .  -  J,  is  explicitly  obtained 

gain  loss 

by  considering  the  dynamics  of  elastic  binary  collisions  and  thus  obtain¬ 
ing  a  differential  cross  section. 


Since  all  the  particles  have  the  same  mass  the  conservation  of 
momentum  law  reads, 

,  *  i  * 

(2-4)  v  4-  v  =  v  +  v 

,  *  T  * 

where  v  and  v  are  the  initial  velocities  and  v  ,  v  are  the  final 

velocities  in  a  binary  collision  of  two  particles. 


We  are  also  considering  only  elastic  collisions  therefore. 


(2-5) 


,,2.1 

ir  1 2 

t  *  1 2  , 

1 1  *.2 

v  + 

1  r>j  •  1 

v  = 

1  rv/  l 

v  + 

1  1 

V 

1  ' 

• 

In  the  lab  frame  the  collision  looks  like 
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By  translating  our  reference  frame  to  a  velocity  V  =  v  4-  '  v  ,  and 


letting 


(2-6) 


A  A 

g  =  v  -  ' v  and  g  =  v  -  1 


V 


we  can  view  the  collision  in  the  center  of  mass  frame,  where  g  is 

* 

the  relative  velocity  of  the  incoming  particles,  and  g  is  the  rela- 
tive  velocity  of  the  outgoing  particles.  Because  it  is  an  elastic 
collision  we  know  that 


(2-7) 


g  =  § 


Pictorically ,  we  have  figure  2,  where  0  ,  and  cf>  are  the  angles  between 
* 

g  and  g  . 


Since  the  motion  of  the  two  particles  before  and  after  the 
collision  is  symmetrical  about  the  center  of  mass  0  ,  we  can  replace 
the  two  particle  scattering  problem  by  an  equivalent  one  particle  scat¬ 
tering  problem,  with  a  fixed  center  located  at  the  center  of  mass.  We 
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use  spherical  coordinates  0  ,  and  <p  ,  with  0  measured  off  the  axis 

through  0  ,  parallel  to  g  ,  where  the  incoming  particle  has  velocity 

■k 

g  and  scatters  off  with  velocity  g  (see  figure  3). 

The  calculation  of  the  collision  term  J  .  -  J,  will  be 

gam  loss 

3 

done  by  considering  an  elemental  volume  d  r  about  the  point  r  ,  and 

focussing  our  attention  on  those  molecules  with  a  velocity  that  is  within 
3 

d  v  about  velocity  v  .  We  will  say  that  a  particle  with  velocity  in 
3 

the  range  d  v  ,  about  v  ,  will  leave  it  if  it  suffers  a  collision  with 

3 

another  particle,  i.e.  we  have  assumed  that  d  v  is  small  enough  that  a 

particle  in  that  range,  suffering  a  collision,  will  leave  it.  In  the 

3  3 

spatial  volume  d  r  there  are  also  molecules  in  the  range  d  fv  ,  about 

some  velocity  ’  v  ,  which  act  as  an  incident  beam  of  particles  on  those 

3 

with  velocities  in  the  range  d  v  ,  about  v  .  The  flux  of  the  incident 
beam  is 


(2-8) 


H(r,'v,t)  d3|v  [ g | 

<*■0 


The  number  of  particles  that  serve  as  "scattering  center"  in  dr  is 
(2-9)  N(r,v,t)  d3v  . 


Define  the  quantity  a(v,'v:v  ,’v  )  such  that  if  "I"  is  the 
incident  flux  then, 


Ia(v,’v:v  ,’v  )dft  =  number  of  particles  deflected  per  sec.  into 

the  solid  angle  dft  (see  figure  3)  centered 
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FIGURE  2. 


FIGURE  3 
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about  the  angle  0  between  g  and  g 

n*/  /'v-' 

JU 

in  the  direction  of  g~ 


a(v,'v:v  ,  fv  )  signifies  collisions  of  the  type  where  particles 


with  velocities  v  ,  and  ’  v  ,  go  into  velocities  v  ,  and  *  v 

3 

The  number  of  collisions  occuring  in  dr  during  the  interval 

*  * 

dt  of  the  type  where  v  ,  and  v  go  into  v  ,  and  ’ v  is 


(2-10) 


The  solid  angle  dfi  is  a  compact  way  of  denoting  the  velocity  ranges 
d  v  and  d  v  about  v  ,  and  ’ v  respectively.  Multiplying 
(2-10)  by  the  number  of  scattering  centers,  (2-9),  we  get 


(2-11) 


This  represents  the  number  of  collisions  that  occur  in  the  time  interval 

3 

dt  of  the  type  where  the  incident  particles  are  in  the  range  d  v  and 


particles  are  in  the  ranges  d  v  and  d  ’  v  centered  about  v  ,  and 
* 

v  respectively.  Since  we  are  interested  only  in  those  particles  with 
velocity  v  that  suffer  collisions,  the  other  incident  particles  'v 
can  have  any  value,  and  the  scattered  particles  can  go  into  any  value. 

Thus  the  number  of  collisions  per  unit  time  occuring  to  particles  in  the 
a3 

range  d  v  is 
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(2-12)  Jlosgd3v  =  d3v  [  d3,v  [  dfi  a (v, ' v: v  ,'v  )|g|  N(r,v,t)  N(r,'v,t)  . 


Equation  2-12  then  tells  us  the  number  of  particles  per  unit  time  that 

3 

leave  our  elemental  volume  d  v  due  to  collisions. 


To  calculate  the  number  of  particles  that  enter  our  elemental 

volume  per  unit  time,  we  perform  similar  arguments  to  the  above,  where 

we  now  seek  collisions  of  the  type  where  the  initial  velocities  can  have 

any  value  and  one  of  the  scattered  particles  can  have  any  velocity  range, 

3 

out  where  the  other  scattered  particle  has  a  range  d  v  about  v  . 

3 

Therefore  the  number  of  particles  that  enter  the  range  d  v  about  v 
per  unit  time  is 


(2-13)  J  .  d3v  =  d3v 
gam 


3  * 

d  'v 


*&»  .L 

A  .  A 


d ft  o (v  ,’v  : v, ’ v) |  g  |  N(r , v  ,t)N(r,'v  ,t) 


From  equations  2-12,  2-13,  it  follows  that 


(2-14) 


(J  .  -  J,  )d  v  = 

gam  loss 


* 


a(v  ,'v  :  v , 1  v)  •  |  g  |  N(r,v  ,t)  x 


Q  J.  A  J. 

x  N(r,’v,t)d  v  d  *  v  d^  - 


"k  k 

a(v, ' v: v  ,'v)  | g|  x 


3  3 

x  N(r,v,t)  N(r,’v,t)d  vd  fv  dft 


We  will  now  assume  microscopic  reversibility  in  collisions  so 


that 
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(2-15) 


A  A  A  A 

o (v  ,  v  :v,  v)  =  a(v,'v:v  ,'v  ) 


Making  use  of  equations  2-7,  2-15,  and  the  fact  that 
3  3.  3*3* 

d  v  d  v  =  d  v  dv  (1)  ,  equation  2-14  reduces  to 


(2-16)  (Jgain  "  Jloss)d3v  =  ^  J  J  ° (v, T  v:  v* , '  v*)  |g|  (N(r,v*,t)  x 


x  n(£>'v  ,t)  -  N(r,’v,t)N(r,v,t)}d3,v  dQ 


Using  equation  2-16,  equation  2-3  finally  yields  the  classical  Boltzmann 
transport  equation,  i.e., 


/9  7\  9N  9N 

(2-17)  — :  +  v  •  — 

at  ~  dr 


+ 


1 

m 


_9N 

9v 


a(v, ’ v: v  ,'v  )  g  (N 


’N^- ’NN)d3'v 


dQ 


where  a(v,'v:v  ,'v  )  corresponds  to  the  differential  cross  section  for 

i  *  * 

a  collision  of  the  type  where  v  ,  and  ?v  go  into  v  ,  and  Tv  ,  and 
N  =  N(r,v,t)  ,  fN  =  N(r,Tv,t)  ,  etc. 


Turning  to  the  Relativistic  case,  we  will  assume  that  these 
are  no  external  forces  except  for  a  self-consistent  background  field  which 
will  take  into  account  gravitation.  To  treat  the  case  of  no  collisions  in 
the  relativistic  case,  consider  an  event  x  ,  and  a  4-momentum  at  x  . 

Also  at  x  consider  some  time-like  unit  vector  n  and  orthogonal  to  it 

y 

can  elemental  3-flat  ds  .  ds  will  serve  as  our  target.  Consider  another 
event  x*  such  that  the  4-vector  from  x  to  x?  is  in  the  direction  of 
P^1  ,  and  at  x*  consider  another  elemental  3-flat  of  the  same  dimensions 
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(size) 

range 


and  orthogonal  to  n  (see  figure  4).  Choosing  some  momentum 

y 

do)  about  on  the  pseudosphere,  we  can  write. 


v  =  N(x,P)  ds  |p^n^|  da) 
v’  =  N (x '  ,P)  ds  |pyn  I  da)  , 

p 


which  implies  that 

(2-18) 


vl  =  N(x* ,P) 
v  N(X,P) 


Therefore  the  ratio,  equation  2-18,  is  independent  of  the  size  of  dco  . 
The  elementary  solid  angle  dco  in  the  direction  P^  is  a  measure  of 
how  much  the  world  lines  fan  out.  For  example  if  dco  =  0  ,  this  would 
imply  that  they  would  not  fan  out  about  P^  ,  (see  figure  5). 


Suppose  we  pick  dco  =  0  ,  then  the  world  lines  will  not  spread, 
and  if  we  don’t  allow  any  collisions  the  world  lines  that  cross  ds  at 
X  will  also  cross  ds  at  x’  >  (figure  4).  This  tells  us  that  at 
and  x’  we  have  the  same  number  of  particles,  i.e. 


So  that 


N(x’ ,P) 

N(x,P) 


1 


■ 
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FIGURE  4 


FIGURE  5 
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But  since  we  know  that  the  above  ratio  (equation  2-19)  is  independent 
of  our  choice  of  dco  then  we  can  write 


(2-20)  N(X,P)  =  N(X',P) 

for  all  choices  of  "dm"  .  If  we  now  let  x*  ->  x  along  the  direction 
of  ,  then  the  condition  equation  2-20  tells  us  that  the  directional 

derivative  of  N  in  the  direction  of  P^  is  zero,  i.e. 

(2-21)  PU  — -  =  0  . 

3  x^ 


Equation  2-21  tells  us  that  if  an  observer  moves  along  with  the  particles, 
he  sees  no  change  in  the  density  function.  We  see  that  equation  2-21  is 
equivalent  to  the  classical  expression  equation  2-2,  if  in  equation  2-2 
the  external  forces  vanish. 


But  just  as  we  saw  earlier,  this  only  works  if  there  are  no 

collisions,  because  otherwise  world  lines  can  leak  out,  and  into  our 

volume  (see  figure  4).  To  account  for  collisions  we  proceed  in  the  same 

fashion  as  for  the  classical  case,  by  replacing  the  right  hand  side  of 

equation  2-21  with  the  phenomenological  collision  tern  J  .  -  J,  , 

1  r  gain  loss 

i.e.  for  collisions  we  have 


3N 


J  . 

gams 


OSS 


(2-22) 


P 


- 
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The  calculation  of  J  .  -  J.  proceeds  in  a  fashion  simi- 

gain  loss 

lar  to  the  way  we  derived,  equation  2-14.  That  is,  we  first  have  the 
conservation  of  4-momentum 


(2-23) 


k 

p  +  'p  =  p  +  ?p 


y 


y  y 


where  we  are  using  normalized  4-momenta,  (see  page  13),  i.e. 

g  P^PV  =  -1  .  The  relative  4-momenta  are 
pv 


(2-24) 


g  =  ’P  -P  ,  g  *  =  TP  *  -  P  ' 
y  y  y  y  y  y 


Note  that  these  are  spacelike  vectors. 


In  a  collision  we  have  the  condition, 


(2-25)  g 


* 


O' 

o 

(“•O 


=  g  =  8 


>  i-e.  g  g 

r* 


y 


*  *y 

Sy  g 


which  follows  simply  from  equations  2-23,  2-24,  and  the  normalization 
condition  P  P^  =  -1  .  If  we  define  the  unit  vector 

y 


(2-26) 


P  5-Ji- 

y 


p  +  ’p 


y 


J  4+f 


,  P 


p  +  ?p 

y  y 


y 


f  *  2 
/  /  ■ 


4+g 


We  see  that  by  equation  2-23 


P  =  P 
y  y 


(2-27) 


■k 
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P  describes  the  motion  of  the  center  of  mass,  i.e.  its  tangent  to  the 

r* 

*5* 

world  line  of  the  centroid  and  its  easily  seen  that  g  and  2  are 

°y  y 

— — 

orthogonal  to  P  .  This  means  that  g  and  g  lie  in  a  3-flat 

U  y  y 

normal  to  P^  ,  in  momentum  space.  The  3-flat  in  question  is  called 
the  centroidal  3-flat,  i.e. 


centroidal  3-flat 


We  translate  to  the  center  of  mass  frame  by  choosing  our  time  axis  to 


be  along  P 

u 


In  the  calculation  of  J  .  -  JT  we  will  define  a  "scatter- 

gam  loss 

ing  amplitude"  which  we  will  relate  to  the  differential  cross  section.  We 


define  the  quantity  in  question,  i.e.  W  ,  by  letting 


(2-28)  W(P,’P:P  ,'P  )  N  (x ,  P )  N(x,’P)  dm  d'u>  dw  d'w  dx  , 


be  the  number  of  binary  collisions  occuring  in  the  4-volume  dx  centered 


about  x^  ,  where  the  initial  particles  are  in  ranges  dto  ,  and  d’co 


about  P^  ,  and  ’  P^  respectively,  and  where  the  final  momenta  are  in 

k  k  k  k 

ranges  dm  ,  and  d'co  ,  about  P^  ,  and  'P  respectively.  One 

should  note  that  this  is  the  relativistic  expression  analogous  to  (2-11) , 
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in  the  classical  case.  To  calculate  J_  ,  we  use  the  same  arguments 

loss 

as  before.  That  is,  we  want  to  know  how  many  particles  in  the  range 
do)  about  P^  ,  leave  due  to  collisions  occuring  in  dr  .  To  do  this 
we  consider  the  number  of  collisions  that  occur  in  dx  with  one  of  the 
initial  particles  in  the  range  dm  about  ,  the  other  initial  parti¬ 

cle  with  arbitrary  value,  and  the  scattered  particles  with  arbitrary 
value.  This  number  is 


(2-29)  N( x,  P)dw  dx 


W(P, fP:P 


* 

*P  ) 


N  (x ,  '  P)  d 1  to  dw  d 


* 


w 


To  calculate  J  .  we  are  interested  in  those  particles 

gam 

which  after  collision,  enter  the  range  dto  about  P  .  We  calculate 

r* 

this  gain  by  considering  the  inverse  collision,  i.e.  we  assume  that  the 
initial  particles  can  have  any  values,  as  well  as  one  of  the  scattered 
particles,  but  that  the  other  scattered  particle  must  go  into  the  range 

&  A 

dw  .  The  "scattering  amplitude"  for  this  is  W(P  ,'P  :P,’P)  .  There¬ 
fore  by  equation  2-28  we  get  that  the  number  of  collisions  satisfying 
the  above  condition  must  be 


(2-30) 


dw  dx 


W(P/,’P  :P,'P)  N(x,P  )  N(x,'P  )  dw  d'w  d'w 


Since  equation  2-22  holds  per  unit  volume  in  phase  space  (2)  ,  (in  our 
case  this  means  "per  unit  volume  per  unit  time  per  unit  volume  in  momen¬ 
tum  space") ,  therefore 
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(2-31)  (J 


gam 


-J-.  )dcodx  =  dcodx 
loss 


"k  "k  k 

W(P  , 'P  :  P, '  P)N 


'  N  dw  d '  oo  d ' 


00 


■  [’ 

-  dcodx 

k 

Making  use  of  the  symmetry  W(P,'P:P 
scopic  reversibility,  equation  2-31, 


W(P, '  P:P"  ,  'P*)  N’N  doT  d'co* 


J.  4 


* 


* 


,fp  )  =  W(p  ,’P  : P , P  )  ,  i.e, 


equation  2-22  finally  becomes 


d'oo  . 


micro- 


(2-32)  Py  — 

8x 

Equation  2-32  is  Relativistic  Boltzmann  Collision  equation  in  the  absence 
of  external  forces  (i.e.  of  a  non-gravitational  nature). 


W(P, ’P:P 


* 

•P  ) 


(N 


vc  k 

'  N  - 


■N’N)  d'oo  doo  d'oo 


Next  we  will  show  that  equation  2-32  is  analogous  to  the 
classical  case,  i.e.  equation  2-17.  We  will  assume  that  the  collision 
process  has  axial  symmetry  about  the  line  of  centers  of  the  colliding 

4. 

particles.  We  define  the  angle  0  between  g  and  g  by  the  relation, 


(2-33) 


Cos  0  - 


h  8 


V 


Since  the  collision  process  in  question  is  elastic  the  condition 
(2-25),  and  (2-27)  must  be  satisfied.  We  also  know  that  just  as  in  the 
classical  case,  an  elastic  collision  is  determined  by  g  and  0  .  Incor¬ 
porating  all  this,  we  can  write 


x 


W ( P , ' P : P  ,'P  )  =  W1(g,0) 


— * 

8  (P 
0JV  y 


-Py)  5(g  -g) 


(2-34) 
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The  scattering  cross  section  cr(g,0)  is  introduced  by  writing 


(2-35)  W(P,'P:PVP  )  = 


(1  +  i  g2) 


2  g  1  a(g,0)  «u(P(1*-P)j)«(g*-g) 


The  reason  for  calling  the  quantity  a(g,0)  the  cross  section 
(3)  is  that,  by  substituting  equation  2-35  into  equation  2-28,  we  get 


(2-36)  a(g,0)N’N  dco  d'w 


tt  i 


~~  A  ^ 

P  )  6(g  -g)  dco  d'co  dx  . 


(1  +  Y)2 


‘k  k 

We  also  know  that  in  an  elastic  collision  satisfying  P  +  'P  =  P  +  ’ P 

y  y  y  y 

we  can  describe  the  initial  momenta  relative  to  some  fixed  tetrad 

(M  ,M„ ,M  ,M. )  in  momentum  space  by  either: 

1  Z  j  a 


(1)  specifying  x*®^  of  P^  and  'x*'®*^  of  *P  (see 


equation  1-18),  as  shown  in  figure  6. 


(2)  specifying  x*®^  °f  P^  ,  and  then,  in  the  centroidal 
3-flat  orthogonal  to  P^  ,  having  some  fixed  orthonormal  triad  from  which 


the  direction  of  g^  is  specified  by  'O  , 
g  itself,  (see  figure  7). 


and  e  ,  as  well  as  of  course. 


A  relationship  among  the  elemental  3-areas  dm  ,  d’co  and  dco 


is  (4) 


12  2  —  2  ft  c\ 

dw  d'd)  =  (1  +  ^  g  )  dco  g  dg  sin  A/  dvde 


(2-37) 
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FIGURE  6 


<3 


FIGURE  7 
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We  can  describe  the  final  momenta  in  exactly  the  same  fashion  as  was 
done  for  the  initial  momenta.  For  the  final  momenta  we  also  have  the 
relation  (4) 


(2-38) 


dm  d'co  =  (1  + 


1 


'2)2 


— * 
do 


dg  sin^l 


.  a- 

-r  d« 


de 


Subsituting  equation  2-38  into  the  expression  2-36  we  have, 


(2-39) 


N'N  a(g,0) 


dto  d '  oo 


g(l  + 


1 

4 


1 

2 


g  a 


4 


* 


')  6(g  -g)  dg  x 


x  6co(P  ~-P  )dw 

y  y 


-n5'  a1  'k 

0  d<v  de  dx 


sm  r. 


It  is  seen  that  the  only  "time" 

—  k  — 

pseudosphere  when  P  =  P  , 

y  y 


(2-38)  is  non- vanishing 

A 

and  when  g  =  g  .  2-39 


is  on  the  unit 
thereby  becomes 


(2-40)  ga(g,0)  N'N  dm  d'to  sin  r{)  d$  de  dx 

Replacing  the  solid  angle  sin  /  of  de  by  d$I  we  can  then 
say  that  the  number  of  collisions  with  relative  3-momenta  g  ,  involving 
particles  with  initial  momenta  in  ranges  dm  ,  and  d'co  ,  that  scatter 
into  the  solid  angle  dQ  is  (5) 


(2-41) 


g  a(g,0)  N'N  dto  d'co  dfi  dx  . 
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We  can  thus  rewrite  the  Boltzmann  collision  equation  to  get, 


(2-42) 


ai  9N 


3x 


U 


g  cr  (g 


&  k 

, 0)  ('N  N 


-  ’NN)  dfi  d’u  . 


Equation  2-42  is  the  relativistic  analog  of  the  classical  Boltzmann 
collision  equation,  (2-17). 
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CHAPTER  III 

Conservation  Laws 


In  this  chapter  certain  general  conservation  laws  will  be 
shown  to  be  direct  consequences  of  the  Boltzmann  Collision  Equation. 

We  will  also  show  how  one  can  obtain  Numerical  flux  vector  and  Energy- 
Momentum  tensor  of  a  Perfect  Fluid  with  explicit  representations  for 
P  ,p  ,  and  y  . 

As  was  mentioned  earlier  a  non-equilibrium  situation  will 
tend  to  equilibrium  by  a  reordering  if  you  will,  of  non-uniformities. 

In  so  doing  certain  quantities  will  be  transported,  e.g.  heat  will  be 
conducted  in  order  to  smooth  out  temperature  differences,  momentum  flows 
in  order  to  smooth  out  pressure  differences,  "mass"  flows  in  order  to 
smooth  out  density  differences,  etc.  These  flows  or  transport  effects 
can  be  described  once  we  have  solved  the  Boltzmann  collision  equation 
for  the  distribution  function  N  .  But  there  are  certain  facts  which  can 
be  deduced  about  these  flows  without  explicitly  solving  the  Boltzmann 
equation.  These  facts  take  the  form  of  conservation  laws.  Let  us  first 
consider  the  classical  case. 

Consider  a  binary  collision  taking  place  at  some  point  r  . 

Let  \p(c.v)  be  any  conserved  quantity  associated  with  a  particle  whose 

v  i .  e .  , 


velocity  is 
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(3-1) 


'k  "k 

ip  +  *  \p  -  ip  +  1  \p 


Consider  the  following  expression 


(3-2) 


JW  = 


* 


a(v,'v:v  ,’v  )|g  ^ ( T N  N  - 


fNN)dfi  d3v  d 3 1 v 


J  (ij))  is  just  the  collision  terra  in  the  Boltzmann  equation  multiplied 
by  i|)  and  integrated  over  v  .  We  now  write  three  equivalent  expres¬ 
sions  of  the  above  integral  by  performing  a  series  of  transformations  of 
the  variables  of  integration,  and  making  use  of  some  symmetries  of  a  . 
Firstly,  we  let  v  ->  ’ v  ,  and  ' v  ->  v  ,  then  equation  3-2  becomes 


(3-3) 


J(»  = 


a('v,v:v  'v  )|g|  '  ip('N  N  -  '  NN)  dft 


j  3 
d  v 


v 


By  means  of  rather  obvious  symmetry  property,  i.e. 


*  * 
a (v, ' v: v  , ’ v  ) 


A 

=  a(’ v,v: v  '  v  ) 


> 


equation  3-3  becomes 


(3-4)  J(» 


a (v* v:v 


v*)  |  g  |  'i/K'n'V 


r  NN)  dft 


,3  j3, 

d  v  d  v 


By  next  performing  the  following  relabeling  of  variables, 

k  k 

v  ->  V  V  ■>  V 

*  * 

' v  ->•  * V  *  v  v 
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from  equation  4-2  we  get 


(3-5) 


JU 

J(f) 


o('v  ,v  :v, 1 v) 


o(v, ' v:v 


i|T(N’N  -  N^'N*)dfi  d3v"  d3,v* 

JU  ^  o  o 

ip  ( ' N  N  -  N'N)  dft  d  v  d  'v  , 


where  X\re  have  made  use  of  microscopic  reversibility,  Liouville’s  theorem 
3  3  3  *  3  * 

(i.e.  dvd*v=dv  d’v)  ,  and  equation  2-7.  By  similar  arguments 
it  can  be  shown  that 


* 


(3-6)  J(>  )  =  - 


r 


^  ^  A  O  O 

o(v,'v:v  ,'v")  'i|f  (? N  N  -N’N)  dfl  d  v  d  'v 


Since  equations  3-2, 4, 5, 6  are  all  equivalent  expressions  then 


clearly 


(3-7) 


J(i|0  =  j  (J(ijj)  +  J(’^)  +  J(’iP  )  +  J(ip  ))  . 


Substituting  equations  3-2, 4, 5, 6  into  equation  3-7  we  get, 


(3-8)  J(<|0  = 


f  ^  vt  JL  O  O 

o(v,’v:v  ,’v  )('N  N  - T  NN)  (i|H-’  ip-ip  -'ip  )  dfi  d  v  d  'v 


But  \p  is  a  collisional  invariant,  so  that  J(ip)  =  0  .  It 
then  follows  from  equation  3-2,  and  equation  2-17,  that 


^ (r~  +  v  •  +  —  F  •  ~ -)  N(r,v,t)  d3v  =  0  . 

3t  ~  3r  m  ~  3v  ~  ~ 


(3-9) 


1 


■ 
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Equation  3-9  contains  the  crux  of  our  conservation  lav/s  and  transport 
effects.  Writing  equation  3-9  out  explicitly  and  then  adding  and 
subtracting  similar  terms  we  have 


inv  n  ,  3N  ,3  ,  .  3N  ,1  ..  3N  .3 

(3-10)  0  =  U  —  d  v  +  U  v.  - - +  —  hpF.  - —  d  v 

J  9t  J  l  3x.  mj  l  3v, 


i  =  1,2,3 


a 

at 


ipN  d  v  + 


.  3N  a3  . 

v.  — -  d  v  + 

1  dX, 

1 


V.  —  N  d3v  - 

1  dX. 

1 


dip  ,3 

v  — - N  d  v 

1  d  X . 

1 


1 

& 

<ro 

.3  1 

• 

3F. 

l  ,3 

1 

+  — 
m  j 

ip  F. 

l  9x . 
i 

d  v  H - 

m 

y 

N  d  v  - 
9x. 

l 

m 

y 

8F. 

ip  N  d  v 

dX. 

1 


+ 


m 


dip  ^  ,3  1  I  „  a^  „  j 3 

F.  — — —  N  d  v - F.  -rm—  N  d  v  . 

l  9v.  m  l  dv. 

l  J  l 


Collecting  terms  equation  3--10  becomes 


(3-11) 


IE  *  N  d3v  +  3*. 


(ipv.N)  d3v  - 


9^  M  ,3  , 

v.  N  d  v  + 


9x .  i 

l 


+  -  f  (4*F .  N)  d3v  -  - 

m  9v.  l  ra 


1 

m 


l 

9F 


^  F.  d3v 


9v.  i 

l 


* 


9v. 

l 


i  M  ,  3 
N  d  v 


0 


By  Gauss's  theorem  the  fourth  term  becomes 


1 

m 


(^F  N)  d3v  =  ^ 
9v.  l  m 


ds  n*  (t^FN) 


-y-yco 


But  because  of  the  fact  that  N  ->  0  as  v  -*  00  we  can  neglect  this  term. 
Defining  the  mean  with  respect  to  the  distribution  function  N  ,  of  a 
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quantity  A  to  be, 

(3-12) 


A  E 


A  N  d  v 


N  d3v 


=  —  A  N  d  v 


It  must  be  pointed  out  that  A  is  really  a  local  mean,  i.e.  its  an 
average  over  velocity  space  but  it  varies  from  point  to  point  in  posi¬ 
tion  space.  Using  equation  3-12,  equation  3-11  becomes 


(3-13) 


9 

3t 


(nip)  + 


9x . 

l 


(mjjv^)  -  n  (v^ 


-ii>  -  s.  (f 

3x.'  m  K  1 


dV.  m  dv. 


4>)  =  o 


Let  us  for  the  moment  consider  the  situation  where  F  is  inde- 

/n; 

pendent  of  velocity,  equation  3-13  then  becomes 


(3-14) 


a 

3t 


(n^) 


+ 


9x . 

l 


(wpv.)  -  n  (v. 
l  l 


i  (F. 
dx .  m  l 


dip 

dv. 

1 


)  =  0 


The  collisional  invariants  we  have  are: 


(1) 

(2) 

(3) 

(4) 

where  we  call 


\p  =  constant,  i.e.  ip  =  1  (for  convenience). 

*  ■k 

conservation  of  mass,  i.e.  m  +  ?m  =  m  +  Tm 

conservation  of  momentum,  i.e.  ip  =  mv^  . 

conservation  of  thermal  energy,  i.e.  ip  =  m  |  v 

v  -  v  the  thermal  velocity. 


The  four  collisional  invariants,  by  means  of  equation  3-14 

give  rise  to  4  conservation  laws. 

for  n  ,  i.e.  •—  +  V*(nv)  =  0  . 

’  9 1  ~  ~ 


For  ip  -  1  we  get  the  continuity  equation 
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For  ip  =  m  it  yields  the  continuity  equation  i.e.  3-14  becomes, 


~  (inn)  +  - — 
d  t  9x . 


(nv^m)  -  n  (v. 


9m 
9x . 

l 


)  - 


l  (F  =  0 

m  l  dv, 


which  reduces  to 


(3-15) 


9  (nm) 


9 1 


9x. 

l 


(nm  v^)  =  0 


Defining  density  of  mass  as  p  =  mn  ,  equation  3-15  becomes 


(3-16) 


|£  +  V-(pv)  =  0 


continuity  equation. 


Putting  ip  =  mv  ,  equation  3-14  becomes 


(3-17) 


9 

9t 


(pvj) 


+ 


9x. 

l 


(pViVj  )  “ 


£ 

m 


(Fi 


6..)  = 

ij 


0 


To  get  equation  3-17  into  a  more  familiar  form  we  can  use  the  fact  that 


(3-18) 


(v.v.)  =  (v.  -  v.)(v.  -  V.)  +  v.v. 
1  J  1  1  J  3  1  3 


Defining  a  Pressure  Tensor 


P.  . 
ij 


to  be 


P.  . 
ij 


p(Vp  -  v_L)(v^.  -  v  ) 


9 


and  making  use  of  equation  3-18,  19,  equation  3-17  becomes 
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(r-r  +  V.  “-)  V.  =  -  F. 
3 1  1  3x.  j  m  j 


P  3xi  ij 


or 


(3-20) 


p  (—  +  V*  V)  v  =  —  F  -  V  •  P  , 
dt  ~  ~  m  ~  ~ 


which  is  just  another  form  of  Newton's  2nd  law. 


1  i  — .  2 

Finally  for  \{j  =  —  m  |v  -  v[  ,  equation  3-14  becomes 


1  3  /  i  — 1 2  ,  3  ,  |  — 1 2  1  3  i  — 1 2, 

(3-21)  —  —  (p  v-v  +- —  (pv.  v-v  -  y  P  (v.  - —  v-v  )  =  0 

2  3t  3x.  i'~  2  l  3x.  J 

l  l 


Defining  temperature  by, 


(3-22) 


KT  E  0  5  -  m  (  |  v—  v  I  2 )  , 


and  heat  flux  by 


(3-23) 


1  —  i  — i  2 

q  =  —  m  p  ((v-v) | v-v |  )  , 

r>*>  Z. 


equation  3-21  reduces  to 


(3-24) 


ITt  (pe)  +SL  +  i  air  (p0vi)  +  m  pi:  5  "  0  > 


where  we  have  used  that  fact  that  p  (v_^(v  — v^))  -  P_^_.  .  Noting  that 

P  is  symmetric,  equation  3-24  can  be  written  as. 
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(3-25) 


3  _3_ 
2  3t 


Defining  A..  =  —  m  (r— ^  + 
ij  2  3x . 


■^p-)  ,  equation  3-25  becomes 


j 


i 


(3-26) 


At  this  point  we  should  note  that  the  preceeding  ideas  of  how 
to  relate  heat  flow,  temperature,  etc.  to  microscopic  variables  must 
be  compared  with  experiment.  But  the  comparison  with  experiment  can’t 
be  made  until  we  have  some  notion  as  to  what  the  distribution  function 
should  be. 

The  calculation  of  the  conservation  laws  becomes  much  simpler 
in  the  relativistic  case.  We  start  by  defining 


(3-27)  J(ip) 


where  ^(x,P)  is  any  tensor  function  of  x^  ,  and  .  Then  in 


exactly  the  same  fashion  as  equation  3-8  was  derived,  we  can  prove  that 


(3-28)  J(ip)=j 


In  the  relativistic  case  there  are  only  two  collision  invari¬ 


ants  : 


■ 
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(3-29) 


if  =  i 


^  =  p 


y 


The  apparent  simplicity  in  the  number  of  collision  invariants  is  essen- 
tailly  due  to  our  space-time  approach.  In  the  classical  case  we  had 
momentum,  and  "thermal  energy"  as  collisional  invariants.  But  now  the. 
4-momentum  vector  P^  incorporates  both  these  conditions.  The  other 
collisional  invariant  of  mass  is  taken  into  account  by  the  restriction 
that  all  the  4-momenta  vectors  must  terminate  on  the  same  fixed  pseudo¬ 
sphere  . 

From  equation  3-28  we  see  that  when  ip  =  1  or  P^  we  have, 

(3-30)  J(ip)  =  0  . 


Thus  in  view  of  equation  3-30  we  see  that  equation  2-32  leads  to  the 
result 


(3-31) 

Equation  3-31  is  the 
For  ip  =  1 


\p  PV  du  =  0 
BxV 


relativistic  analog  of  equation 


equation  3-31  yields. 


Dv  BN  ,  n 
P  - -  do)  =  0  , 


3x 


v 


3-9. 


or, 


i  'i 
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(3-32) 


3 


PV  N 


dw  =  0 


For 


equation  3-31  yields 


3N 


3x 


v 


dw 


0  , 


or, 


(3-33) 


P|J  PV  N  dw  =  0 


Noting  that  the  Energy-Momentum  tensor  is, 


T 


pv 


2 

m  c 


PU  PV  N  dw 


) 


and  the  Numerical  flux  vector  is. 


My  = 


m  c 


PP  N  dw 


we  see  that  in  general  equations  3-32,  33  become 


(3-34) 


0 


(3-35) 


Equations  3-34,  35  become  our  conservation  laws. 
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Perfect  Fluid 


As  an  application  of  the  preceeding  ideas  the  Numerical  flux 
vector,  and  the  Energy-Momentum  tensor  of  a  perfect  fluid  will  be 
derived,  with  the  explicit  expressions  for  the  pressure,  energy  density, 
and  density  of  proper  mass. 

Before  we  do  this  it  will  be  necessary  to  derive  the  equili¬ 
brium  distribution  function.  The  relativistic  derivation  of  the  equili¬ 
brium  distribution  is  similar  to  that  of  the  classical  case.  That  is, 
we  first  formulate  the  relativistic  H-theorem  and  then  use  the  equilibrium 
condition  to  obtain  the  equilibrium  distribution. 

We  define  the  entropy- flux  vector  to  be 


(3-36) 


Sy  =  -  K 


N(x,P)  log  N(x,P)  Py  dm 


where  K  is  the  Boltzmann  constant.  Suppose  that  in  equation  3-28  we 
put  ip  =  log  N  ,  then  we  get 


(3-37) 


JOlO  =  J  K 


ju  at  ^  ^ 

W(P,'P:P  ,?P  ) (log  N  +  log  TN  -  log  n"  -  log  'N  ) 


^  ^  A  ^ 

(TN  N  -  ,NN)d'm  doo  d'mdm 


K 


W(P,  TP:P*,  rP~)  (log(’N/'N*)  -  log ( * NN) ) 


V.U  JU  JU  ^ 

(fN  N  -  ,NN)d,co  dm  d’w  dm  . 
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But  we  also  know  that  from  equation  2-31,  27,  we  have 


(3-38) 


log  N  Py  doo  =  J(i/>) 


3x 


y 


It  can  also  be  seen  that 


(3-39)  Sy  |  - - —  N  log  N  Py  doo  =  -  log  N  Py  doo 

IP  J  J 


3x 


9x 


_y  3N  , 

P  -  doo 


3x 


y 


By  equation  3-32,  we  see  that  equation  3-39  simp lies  to 


(3-40) 


,y 


y 


pP  log  N 


3N 


3x 


y 


doo 


Noting  equations  3.37,  and  3-38,  equation  3-40  becomes 


(3-41)  S 


y 


1  rrr 


y  4 


K 


W(P  ,  'P  :P~,  TP~)  [log ( 'N~N  )  -  log ( ’ NN) ]  x 


x  [’N  N  -’NNjd’oo''  doo  d'oo  doo  . 


A 

In  equation  3-41  consider  the  expression  [log(rN  N  )  -  log(’NN)]  x 


* 


[ 'N  N  -’NN]  ,  where  we  know  that  N,’N,N  ,  and  'N  2l  0  .  This  expres¬ 


sion  is  of  the  form  (x-y)  log  (— )  .  Consider  the  following  3  cases: 


(1)  if  x  >  y  then,  (x-y)  >  0  a  log  (^-)  >  0  ,  therefore 


(x-y)  log  >  0  . 


(2)  if  x  =  y  then  (x-y)  log  (^-)  =  0 
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(3)  if  x  <  y  then  (x-y)  <0  a  log  (— )  <  0  ,  thus  (x-y) 


log  (-)  >  0  . 

y 


The  result  of  this  being  that 


(3-42) 


V 


x 


>_  0  V  y  >_  0  (x-y)  log  0^-)  >_  0 


Using  equation  3-42,  and  noting  that  W(P,’P:P  ,'P  )  is  positive  definite 
we  see  that  equation  3-41  yields 


(3-43)  Sy|  >0 

IP  ~ 

which  is  Boltzmann’s  H-Theorem.  The  connection  between  equations  3-43, 
equations  1-31  and  the  classical  H-Theorem  can  be  seen  by  letting  ny 
be  coincident  to  Py  ,  and  then  picking  the  time  axis  along  Py  (1) . 

We  know  that  the  equilibrium  condition  is  characterized  by 

(3-44)  Sy  |  =  0  . 

lu 

From  equation  3-41  we  can  see  that 

(3-45)  Sy  |^  =  0  iff  log  N  +  log  'N  =  log  n'  +  log  'N 

k  k 

The  equation  log  N  +  log  ’N  =  log  N  4-  log  TN  represents  some  colli- 
sional  invariance.  But  we  know  that  the  only  collisional  invariants  are 
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P  and  a  constant.  Therefore  log  N  must  be  some  linear  combination 

of  Py  and  a  constant.  Letting  the  constant  be  log  (x)  \>7e  have 


(3-46) 


log  Nq(x,P)  =  loga(x)  +  3^(x)  P 


where  Nq  denotes  the  equilibrium  distribution.  Therefore  we  can 
write 


3  (x)  P1 


(3-47) 


Nq(x,P)  =  a(x)  e 


where  3^(x)  must  be  timelike  so  that  3^(x)PM  _<  0  ,  and  thus  guaran¬ 


teeing  that 


N  dm  is  finite.  Letting 
o 


32  =  -3  (x)3y(x)  ,  and  3  (x) =  7  3  u  , 

y  y  y 


then  u  uy  =  -c2  ,  and  by  substituting  3-47  into  the  collision  term 
P 

of  the  Botzmann  equation  and  noting  the  collisional  invariance  of  Py  , 
we  have 


(3-48) 


8N 

Py  — 2-=  0 


8x 


P 


Substituting  equation  3-47  into  equation  3-48  we  get. 


3  PV  ~  3  PV 

pU  (a(x) e  V  )  =  Py  ( — ^  +  a(x)3v|^PV)e  V 


=  0 


From  which  we  can  write 


Py  +  a(x)  PyPV  3  |  =  0  . 

3xy  Vly 


Since  P  can  be  any  timelike  unit  vector  we  have 


(3-49)  a(x)  =  const.  and  3  I  +  3  |  -  0 

v  |  y  y  |  v 


Thus  equation  3-49  indicates  that  equilibrium  implies  the  existence 


of  a  time-like  Killing  vector  field. 


We  are  now  ready  to  derive  the  required  expression  for  the 
perfect  fluid.  By  defining  the  function 


(3-50) 


3  P 


y 


Z(a,3)=ae  dm  , 
o  y  I 


we  see  that 


(3-51) 


az 


o 


33 


Py  a  e  y  dm  =  Py  Nq  dr 


and 


(3-52) 


32Z 


o 


33  93 

y  v 


3  ph 

PyPV  a  e  y  dm  = 


PyPV  N  dm 
o 


Noting  equations  3-34,  35,  equations  3-51,  52  yield 
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(3-53) 


9Z  .  32Z 

„  y  o  m  yv  2  o 

M  —  me  -  T  —  me - 

0  93  >  Lo  93  93 

P  y  v 


Picking  polar  coordinates  with  8^  parallel  to  the  axis  off  which  x 
is  measured,  equation  3-50  becomes 


(3-54)  Z  (a,  8  )  =  a 
o  y 


■°°  fir  r2  7i  -8  P^1  _ 

y  .  .  2 


e  sinh  x  dx  sin0  d0  d<J> 


o  o 


o 


=  a 


r 

•IT 

V - 

o 

O 

2  7T 


O 


-8  cosh  x  •  r2  ,  -  ,, 

e  sinh  x  dx  sin  0  d0  dcj> 


Using  the  formula 


8 


n 


(3-55)  K  (8)  =  —  .. 

n  1 *3 • • • (2n-l) 


fCQ 

-8  cosh  x  •  i  2n  . 
e  sinh  x  dx 

o 


where  K  (8)  is  the  Modified  Eessel  function  of  the  second  order; 
n 

it  is  seen  that 


Kx(8) 

“ 8  ' 


-8  cosh  x  •  u2  j 
e  sinh  x  dx 


Integrating  over  0  ,  and  <p  ,  equation  3-54  becomes 


(3-56) 


K  (3) 

Va’V  =  ^  “  ~r~ 


From  equation  3-53  we  have 


a  K  (3) 

w  y  9  f ,  i  \ 

Mo  =  ra  c  (4lT  01 
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noting  that 


8  3  3  3 


33  33  33 

y  y 


1  u  3 
c  m  33 


the  above  expression  becomes 


,,  P  ,  1  P,  3  ,Kll6\ 

Mo  =4.a.c(--Oi  (— b~ ) 


By  means  of  the  identity 


ie  [6"n  Kn(6)]  =  -B-n  Kn+1(S) 


equation  3-53  finally  becomes 


(3-57) 


K2(3) 

M  =  4tt  a  m  - -r—  u^ 

o  3 


k2(3) 

Letting  the  density  p  =  4tt  a  m  3  we  see  that  equation  3-57  is 
the  usual  expression  for  the  numberical  flux  vector  of  a  perfect  fluid, 

i.e. 


(3-58) 


y  y 

M  M  =  p  u 
o 


In  a  similar  fashion  starting  with, 


2  2  d 

m  yv  l  o  .  3  3  (— - — ) 

h  =  m  c  =  4rr  a  m  c  e 


y  v 


y  v 
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we  get 


yv  _  ,  K3(3)  yv  2  K2(3)  yv 

(j-jy'  -  47T  ra  a  — - —  u  u  +  4tt  m  a  c  — - — -  g 


3 


3 


where  x^e've  used  the  fact  that  =  g^V 


33. 


Defining  yc  to  be 


v 


the  energy  density,  and  P  for  pressure  and  letting 


p  K  (3)  K  (3) 

(3-60)  y  +  — £  =  4 tt  m  a  — - —  ,  and  — y  =  4 tt  m  a  - y—  , 

c  c  3 


equation  3-59  becomes 


(3-61) 


tq  -  (y  +  — pu  u  +  p  g  , 


which  is  the  Energy-Momentum  tensor  for  a  perfect  fluid.  Defining 

,  .  .  .  y  v 

the  projection  operator  A  to  be 


yv  yv  1  y  v 
A  =  g  +  — ^  u  u  , 


(it  projects  objects  down  onto  the  hypersurface  orthogonal  to  the 
vector  u^)  ,  equation  3-61  can  be  rewritten  as 


(3-62) 


yv  y  v  „  yv 

T  M  =  y  u  u  +  P  Am  . 
o 


It  can  also  be  seen  from  all  this  that  u^  =  is  the  hydrodynamic 

3 

4-velocity.  Relativistic  temperature  is  defined  by 
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(3-63) 


3 


_  m  c 

~  kTt 


From  the  conservation  law  (3-34)  ,  we  get  for  the  perfect 


fluid 


(3-63) 


M  y  =  (puy) ,  =  uP  —  +  P  uy  |  =  0  . 

°  I  11  3X11  ^ 


Defining  the  derivative  along  a  stream  line  -j-^  to  be 


d  dx  3  y  3 

_  =  - - =  u  - 

dx  dx  y  „  y 

3x  3x 


we  finally  arrive  at  the  continuity  equation 


(3-64) 


-  0 

dx  |  y 


The  conservation  law  T  XV]  =0  becomes 

o  I v 


T  XV  -  =  0  =  A  T 

o  v  Ay  o  |v 


But 


T 


Av 


o 


v 


(3^y)u\iV  + 


y  u 


v 


V 


u 


+  y 


A  v 
u  u 


+  (3  P) A 
v  v 


Av 


+  P  A 


Av 


v 


+  (3  P)  AXv 
V  7 


9 


62 


Therefore  we  have 


or 


A  T  XV, 

A  y  o  v 


p  6u  , 

(y  +  -t)  +  A  3,  p  =  o 

2  6t  y  A 

c 


(3-65) 


p  5u  i 

(y  +  ^  =  -  AX  3  P 

c2  fix  y  A 


Equation  3-65  is  essentially  the  relativistic  analog  of 


Newtons  second  law,  where  (y  H - is  the  effective  mass. 
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CHAPTER  IV 

Non-Equilibrium  and  Transport  Effects 


In  this  chapter  we  investigate  transport  effects  arising  from 
non-equilibrium.  By  means  of  a  relativistic  version  of  the  Grad  Method 
we  obtain  explicit  representations  of  the  Transport  coefficients. 


In  the  last  chapter  it  was  shown  that  for  the  equilibrium 
situation  we  obtained  the  Energy-Moment urn  tensor,  and  the  Numerical  flux 
vector  of  a  perfect  fluid.  In  so  doing  we  identified  u^(=c6^/g)  with 
the  hydrodynamic  4-velocity.  This  identification  was  unique  in  that 
u^  was  naturally  singled  out  as  a  result  of  the  equilibrium  condition. 

As  a  consequence  the  Energy-Momentum  tensor,  and  the  Numerical  Flux 
vector  we  expressed  entirely  in  terms  of  u^  ,  and  certain  scalars. 

But  given  a  non-equilibrium  situation  the  choice  of  the  Hydrodynamic  4- 
velocity  becomes  somewhat  arbitrary  (1).  Because  once  we’ve  lost  equili¬ 
brium,  we've  lost  the  condition  that  singled  out  (3^  .  Since  we  will  be 
interested  in  viscous  effects  and  the  associated  transport  coefficients 
it  will  be  necessary  to  retain  conventional  Hydrodynamics.  This  is  done 
by  assuming  that  the  deviation  from  equilibrium  of  the  system  is  small, 
so  that  we  can  say  that  at  each  point  in  space-time  the  distribution 
function  can  be  approximated  by  a  local  Maxwellian  distribution, 


Nq  (x,P) 


a 


(x) 


Vx)p 


(2) 


« 
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From  this  it  can  be  seen  that  No(x,P)  will  vary  from  point  to  point, 
which  means  that  u^  will  also  be  a  function  of  x  . 

The  problem  still  remains  as  to  how  Nq(x,P)  will  be  fitted 
to  the  actual  distribution  function.  This  will  involve  removing  a  certain 
amount  of  the  arbitrariness  in  the  five  parameters  a(x)  ,  B^(x)  occuring 
in  Nq  .  More  on  this  will  be  said  at  a  later  point  in  the  text.  The 
condition  of  small  deviations  from  equilibrium  is  expressed  as, 

(4-1)  N  (x , P)  =  Nq(1  +  f  (x , P )  )  , 

where  f(x,P)  (<<1)  is  the  first  order  correction  to  the  equilibrium 
distribution.  We  will  therefore  neglect  terms  second  order  and  up  in  f  . 
We  also  impose  the  condition  that 

3  f  3  N 

v  «  _v_ 

f  N  * 

Substituting  equation  4-1  into  the  expression  for  the  energy-momentum 
tensor  we  get 


uv  2 
T  =  me 


N  (l+f)PyPV  dm 
o 


=  me 


r\ 

N  PyPV  da)  +  me 
o 


N  f  PyPV  do)  , 
o 


which  leads  to 


TPV  =  T  uv  +  T  yv 
o  1 


J 


(4-2) 
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where 


(4-3) 


T 


2 


N  f  PyPV  dw 


1 


me 


o 


Similarly  for  the  Numerical  flux  vector  we  get. 


(4-4) 


My  =  M  y  +  M.U  , 
o  1 


where 


(4-5) 


M, V  =  me  N  f  PP  do) 


o 


and  T^yV  will  contain 


the  necessary  transport  effects.  To  get  a  better  understanding  as  to  the 
physics  involved,  in  the  preceding  statement,  let  us  for  the  moment 
approach  of  non-equilibrium  from  a  more  phenomenological  point  of  view. 

The  ideas  we  will  use  were  first  put  forth  by  Eckart  (3) ,  and  more  recent¬ 
ly  elaborated  on  by  Weinberg  (4).  In  the  Eckart  approach  the  arbitrariness 
involved  in  the  choice  of  the  Hydrodynamic  4-velocity  is  removed  by  stipu¬ 
lating  that  the  density  p  is  the  4th  component  of  My  in  the  comoving 
frame,  in  other  words 


(4-4) 


y 


also  that  the  energy  density  y  is  the  (4,4)  component  of  T 
measured  in  the  comoving  frame,  in  other  words 


as 
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(4-5) 


u  u  TyV  =  u  u  T  yV  , 
y  v  y  v  o 


and  by  finally  choosing  the  space-time  direction  of  to  be  that  of 

the  particle  current,  i.e.. 


(4-6) 


U 

u  = 


(-lA  ) 

a 


°  My 


1/2 


By  means  of  A^V  ,  and  u^  we  now  proceed  to  decompose  the 
energy-momentum  tensor  into  its  proper  components,  i.e. 


(4-7) 


y  -  u  u  t 
a  p 


a3 


(4-8) 


Q11  5  -Ay  T°B  u, 
a  i 


(4-9) 


At  t“B 

a  8 


From  equations  4-7, 8, 9  we  get  the  following  identity 


(4-10) 


yv  yv  yv  vy  yv 

T  =yuu  +  Qu  +  Qu  +  e 


The  components  of  (equations  4-7, 8, 9)  have  an  important 

physical  significance,  and  through  them  equation  4-10  reveals  the  struc¬ 
ture  of  T^V  .  To  make  this  physical  significance  more  transparent  let 
us  confine  our  analysis  to  the  local  rest  frame  i.e.  choose  the  time 
axis  to  be  along  u^  .  Then  using  equation  4-5  we  have  that  y  as 
defined  by  equation  4-7  becomes 


■ 
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(4-11) 


y  = 


44 

T 

o 


which  is  just  the  invariant  energy  density. 

Defining  what  is  commonly  designated  as  the  viscous  stress 
tensor  (5) ,  xyV  ,  through  the  relation 


(4-12) 


yv 

e 


xyV  +  PAyV 


» 


where  P  is  the  local  equilibrium  pressure,  equation  4-10  becomes 


(4-13) 


:yv  =  T  yV  +  QV  +  QV 

o  x 


+  X 


yv 


where  we  have  used  equation  3-61.  Comparing  equation  4-13  to  equation 

4-2  we  now  see  that  the  phenomenological  addition  T^yV  to  the  perfect 

fluid  energy  momentum  tensor  T  yV  is 

o 


(4-14) 


,yV  =  QUuV  +  QVuy  +  xyV 


which  includes  heat  conduction,  and  dissipative  effects  due  to  viscosity. 

To  get  more  explicit  representations  of  Qy  ,  and  xyV  let  us  once 
again  confine  our  analysis  to  the  local  rest  frame.  The  usual  phenomeno¬ 
logical  arguments  given  in  ascertaining  the  general  form  of  x1'*  (i,j  ,=1,2,3) 
are  essentially  (5,6)  : 

(a)  x1''  must  depend  on  the  space  derivatives  of  velocity,  since 

its  the  relative  motion  of  the  fluid  parts  that  cause  inter¬ 
nal  friction  ,  i.e.  viscosity. 


68  - 


(b)  The  vel  ocity  gradients  are  assumed  small,  so  that  we  can 
expect  that  x1^  depends  on  the  first  derivatives. 


(c)  cannot  be  independent  of 


8  u 

9x* 


since  x1'*  =  0  when 


u  =  const. 


(d)  As  a  consequence  of  b  ,  the  derivatives  should  appear 
linearly. 

(e)  If  the  fluid  is  in  uniform  rotation  then  there  will  be  no 
relative  motion  hence  no  viscosity.  We  know  that  for 
uniform  rotation  that  the  rate-of-strain  tensor 


3u 


r  + 


3  uJ 


=  0  .  Thus  sums  of  the  form 


3  u  .  3  u~ 


—  + 


should 


3xJ  3x  3xJ  3x 

appear.  Note  that  this  expression  satisfies  conditions 

a  ,b  ,  c  ,d  . 

Hence  the  most  general  form  of  x1^*  satisfying  conditions  a,-d 


is 


(4-15) 


T 


ij  .3u  3uJ 

=  -  v  ( - r  +  - ; 


V-u  d1^)  -  xV-ud1^ 


3xJ 


3x 


/v;  /vi 


v  is  the  shear  viscosity,  and  <  is  the  bulk  viscosity.  The  extra 

term  of  -  \  V  •  u  d^  in  equation  4-15  is  put  in  just  to  make  the  quan- 
J  ~  ~ 


tity  in  paranthesis  traceless.  In  general 


(4-16) 


uv  _ 


.p  . v  ,  a  3  ,  3  a  2  a3  A  >. 

-  v  A  A  0  (u  1  +  u  1  -  t  A  u  >  ) 
a  3  d  A 


mv  X 
K  A  u 
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For  future  reference  we  make  the  following  definition 


(4-17) 


0PV  =  AVa(u°lB  +  uBla 

a  3 


2  .  ap  X  N 

A  u  .  ) 


Equation  4-16  now  becomes 


(4-18) 


yv  yv  Ayv  X 

x  =  -v  a  -  k  A  u 


X  * 


To  obtain  the  phenomenological  form  of  the  heat  flux  vector 

we  once  again  refer  to  the  arguments  given  by  Weinberg  (  )  .  To  obtain 

the  structure  of  we  once  again  for  simplicity  consider  the  problem 

in  the  locally  co-moving  Lorentz  frame  i.e.  u^  =  (0,0, 0,c)  ,  Thus 

i  i4 

from  equation  (4-14)  we  can  see  that  Q  (i  =  1,2,3)  is  just  T^ 

The  structure  we  pick  for  Q1  must  be  such  that  it  reduces  to  the  classi¬ 
cal  Fourier  conduction  Law  under  suitable  conditions.  The  most  general 

GL 

form  which  will  yield  positive  entropy  generation  (  ),  i.e.  S  | ^  >_  0  , 
is 


Q1  =  T^4 


Thus  in  general  we  can  write 


(4-19) 


Qy  =  ..xiua  tl3L  +  T  i  j  . 

dX 


Notice  that  for  small  accelerations  4-14  reduces  to  C>  =  -XVT  .  The 
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relativistic  version  of  the  heat  flow  vector  gives  rise  to  a  result  not 

occuring  in  the  classical  case,  i.e.  the  terra  T  u  gives  rise  to 

a 

isothermal  heat  flow  in  accelerated  matter,  in  the  opposite  direction 
to  the  acceleration. 

For  our  future  calculations  we  will  use  an  alternative  (  ) 

definition  of  Qy  given  by  Israel  (  )  : 


(4-20) 


Qy  =  X  Ayv  9  0 
v 


where 


e 


TS)  . 


In  the  preceding  discussions  on  the  phenomenological  structure 
u  yv 

of  Q  and  x  we  followed  the  approach  of  Eckart.  In  so  doing  the 
arbitrariness  in  the  hydrodynamic  4-velocity  was  removed  by  imposing 
equation  4-4, 5, 6.  In  what  will  follow  we  will  not  impose  all  of  these 
conditions.  We  will  still  retain  equations  4-4,5  but  not  equation  4-6, 
i.e.  we  will  not  assume  that  u^  and  have  the  same  space-time  dir¬ 

ection.  As  a  consequence  of  this  the  definition  of  Q  y  will  have  to  be 
modified.  We  will  define  heat  flow  to  be  the  flow  of  energy  relative  to 
the  particle  stream  (8),  i.e. 

Qh  E  -  c  - - - pyx  [TyVM - MyTXv  MM] 

(-M  Ma)1/  Z  v  (-M  M  )  X  v 

a  /  3 


(4-21) 
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Comparing  equation  4-21  with  equation  4-8  we  see  that  equation 
4-21  would  reduce  to  equation  4-8  if  we  were  to  retain  equation  4-6. 
Since  we  will  be  concerned  with  quantities  expressed  only  up  to  first 
order  in  f,Qy  expressed  to  first  order  in  f  is, 

(4-22)  Qy  =  -  n  c2  My  -  Ay  T.pV  u  (9) 

1  p  1  v 

where 

c 

n  =  - 

p 


The  viscous  stress  tensor  will  be  given  by 


(4-23) 


yv 

T 


Ly 


T  aS 

8  1 


Equation  4-23  follows  directly  from  equations  4-9,12. 

The  viscous  and  heat  flow  effects  are  contained  in  M^y  and 
T^yV  .  But  to  get  expressions  for  these  we  still  have  to  solve  for  the 
function  f  .  To  solve  for  f  we  will  consider  the  linearized  Boltz¬ 
mann  equation.  To  obtain  the  linearized  Boltzman  equation  one  substitutes 
into  the  Boltzmann  equation  (i.e.  equation  2-32)  the  first  order  stipula¬ 
tions  associated  with  equation  4-1.  Upon  performing  this  substitution  we 


have , 
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(4-24)  Py9  (N  +  N  f)  = 
y  o  o 


W  { (N  *  +  N  ~f*)('N  *  +  'N  *’f*) 
o  o  o  o 


^  J 

-  (N  +  N  f)('N  +  'N  'f)>  d'w  dm  d'w' 

o  o  o  o 


Expanding  out  equation  4-24  we  have, 


(4-25)  Py 9  N  +  PmN  9  f  +  PMf9  N  = 
y  o  o  y  y  o  j 


>y. 


'k  -k  k  k 

W  (N  'N  +  N  ' N  (f  +  'f 
o  o  o  o 


+  f 2 (N  "'N  k)  -  N  TN  -  N  *N  (f  +  ' f)  -  f 2 (N  fN  )} 


o  o  o  o  o  o 


o  o 


k  k 

x  d  ’  oj  dco  d  '  w 


Since  we  are  only  retaining  terms  up  to  first  order  in  f  equation  4 
reduces  to 


(4-26)  Py9  N  = 

y  o 


A  A  A 

W  {  (N  ' N  -  N  '  N  )  +  N  ' N  (f  +  '  f  )  - 

o  o  o  o  o  o 


N  TN  (f+' f ) }d ' w  dm  d? 
o  o 


But  we  know  that  N  is  the  local  equilibrium  distribution  function 


o 


hence  it  must  satisfy  equation  3-45,  so  that, 


N  k  ’N  k  =  N  'N 
o  o  o  o 


In  view  of  this  equation  4-26  becomes 


(4-27)  Py9  N  = 
y  o 


=  -  N 


o 


W  TN  (f+'f-f*-'f*)  d'w  dm*  d’w* 
o 


* 

)  + 


X 


-25 


k 

w 
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Defining  6(f)  =  f  +  'f  -  f  -  ’f  ,  and  defining  the  linear  integral 
operator  L  to  be 


L(f) 


W(P,'P:P*,fP*)  'N  6(f)  d'o)  dw*  d’w*  , 

o 


equation  4-27  simple  becomes 


(4.28) 


PVI  ^  Nq(x,P)  =  -Nq(x,P)  L(f)  . 


It  can  easily  be  seen  that, 


(4-29) 


No(x,P)G(x,P)L(F)dw 


1  f 
^  JJ 


W 


N  'N 
o  o 


* 


6(F)6(G)  doo  d'w  dco  d* 


w 


Thus  multiplying  equation  4-28  by  "l"  and  integrating  over  doo  we  can 
see  that  inview  of  equation  4-29  we  get, 


(4-30) 


0  , 


where  we  have  used  the  fact  that  6(1)  =  0  .  By  multiplying  equation 
4-28  by  PV  and  integrating  over  dm  we  get 


(4-31) 


) 


where  we’ve  used  equation  4-29,  and  the  fact  that 


6[PV]  =  0  . 
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By  equation  4-28  it  is  seen  that  in  order  to  solve  for 

f(x,P)  we  have  to  specify  the  space-time  gradients  of  a, 8  i.e. 

y 

3^  .  But  in  view  of  equations  4-30,  31  we  see  that  all  these  gradi¬ 

ents  cannot  be  freely  assigned,  i.e.  they  are  subject  to  the  five  condi¬ 
tions  expressed  in  equations  4-30,  31.  The  ramifications  of  equations 
4-30,  31  have  already  been  derived  in  chapter  III,  (see  equations  3-64, 
65). 


To  obtain  a  solution  for  f  we  make  use  of  the  Relativistic 
version  of  the  Grad  approximation  method  (1)  as  was  suggested  by  Anderson 
(11).  What  one  essentially  does  is  to  construct  an  orthogonal  set  of 
polynomials  in  momentum  by  means  of  the  Gram-Schmidt  Orthogonalization 
process,  where  we  use  the  local  equilibrium  distribution  function,  Nq  , 
as  the  weighting  function.  For  example  designating  the  polynomials  by 

ml  l 

H  ,  where  m  is  the  order  of  the  polynomial  and  (£..  ,...,£  ) 

*  r  j  1  m 

is  the  appropriate  tensor  index,  then  we  have 


o 


H  = 


1 


HP  =  Py 


o 

aU  H 


,  HyV  =  pV 


ayVp 


) 


etc. ,  and  where  ay  ,  ayVp  ,  aPV  ,  etc.  are  determined  by  the  orthogon¬ 
ality  conditions  (12).  Once  the  orthogonal  set  is  constructed  we  expand 
N  as 


N  = 


m£. 


A 


(x)  H 


m 


(x,P) 


75 


The  significant  and  crucial  step  of  the  Grad  approximation 

m 

method  is  to  set  all  coefficients,  a p  «  ,  past  a  certain  order 

1  ra 

equal  to  zero.  For  example  suppose  we  set  all  the  coefficients 
m 

p  equal  to  zero  for  m  >  2  ,  then  equation  4-40  can  be  essen- 

1  m 

tially  written  as, 

(4-33)  N  =  N  (A(x)  +  B  (x)  Py) 

o  y 


But  his  just  reproduces  the  perfect  fluid  theory.  Thus  we  are  led  to 

m 

the  approximation  of  setting  all  coefficients,  a p  „  ,  equal  to 

zero  for  m  >_  3  .  The  obvious  question  to  ask  at  this  point  is;  "Why 
this  order  of  approximation  and  not  higher?".  Andersons  answers  this  by 
saying  that,  "the  best  that  one  can  do  at  present  is  to  increase  the 
order  of  the  approximation  and  calculate  the  effect  on  whatever  quanti¬ 
ties  one  is  interested  in  computing.  This  at  least  in  the  Grad  Method  is 
doable  although  the  effort  involved  is  rather  enormous ♦  Alternately,  one 
can  try  to  compare  predictions  of  the  approximation  with  experiment.  At 
least  in  the  classical  case  these  predictions  agree  remarkably  well  with 
experiment".  (13)  Using  this  approximation  and  comparing  equation  4-1 
with  equation  4-32  we  can  write  without  loss  of  generality  that 


f  =  C  (x)  +  C  (x)  Py  +  b  PyPV 

5  y  yv 


(4-34) 


. 
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To  obtain  the  third  moment  equation  which  will  relate  the 
transport  effects  we  first  multiply  equation  4-27  by  PXpy  and  integrate 
over  the  momentum  i.e. 


(4-35)  PXPypV8  N.  =  - 

v  o 


W  ’N  N  6(f)  PXpy  dm  d’m  dm*  d'm* 
o  o 


k  k 

Using  the  symmetries  of  W  and  the  fact  that  N  ’N  =  N  '  N 


o  o  o  o 


equation  4-35  can  be  written  as, 


(4-36)  8 


N  PXPyPV  dm 
o 


1 


W  'N  N  6(f)  6(P  Py)dm  d'm  dm"  d'm 
o  o 


Substituting  equation  4-34  into  equation  4-36  we  get 


(4.37)  3  U 

vo 


Xpv  1  3 

=  — r-  m  c 
4 


W  N  'N  6 (C  (x))  6 (P  Py)  dm  d'm  dm  d'm 
o  o  5 


1  3  f  f 

-  4  m  c  jj 


W  N  'N  C  (x)  6  (PV)  6(P  Py)  dm  d'm  dm  d'm 

O  O  V 


1  ,  f  v  3 
4  baB(x)  OC 


W  N  'N  6(  PaP6)  {(ph11)  du  d'u 

°  °  k  k 

dm  d'm 


where  we’ve  made  use  of  the  following  definition: 


(4-38) 


n  Afiv  _  3 

Li  =  m  c 


PXPyPV  N  dm 


o 


Since  C^(x)  and  PV  are  collisional  invariants  then  equation 


4-37  simplifies  to 
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(4-39) 


3  U  X,JV  =  -  b  x“SX,J 


vo  a$ 


where 


X 


aBXy  _  1  3 


=  -r  ra  c 
4 


W  N  'N  6(pV)  6  (P  Py)  dto  d’oo  dm  d’to  (14)  . 

o  o 


Equation  (4-39)  contains  all  the  necessary  transport  effects.  From  this 

we  will  proceed  to  derive  expressions  for  the  transport  coefficients. 

But  before  proceeding  into  the  explicit  calculations,  it  is  best  that 

the  essential  overall  plan  be  described  so  that  the  reader  can  appreciate 

the  rather  beautiful  simplicity  of  the  results  when  they  do  finally 

appear.  We  will  express  the  left  hand  side  of  equation  4-39  in  terms  of 

gradients  of  thermodynamic  parameters,  and  of  Hydrodynamic  velocity.  In 

particular  we  will  have  terms  like,  a  n  (see  equation  4-17)  ,  AP  0  , 

r  a3  a  ,u 

etc.  Through  equations  4-22,  23  we  will  be  able  to  express  b^  in 

terms  of  t  ^  ,  and  Q  .  Thus  we  transform  equation  4-39  into  linear 
a$  a 

a  O 

combinations  of  o  ^,A  P0  „  ,  etc.  on  the  left  hand  side  and  x  .  ,  and 

a3  ,3  a3 

Q  on  the  right  hand  side.  By  properly  identifying  corresponding  terms 
a 

on  each  side  we  end  up  with 


?  Ayp  a  0 


P 


But  these  are  just  the  phenomenological  expressions  for  the 


viscous  stress  tensor,  and  heat  flux  vector  we've  derived  (equations  4-18, 
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4-20).  Thus  we  can  identify  the  coefficients  (represented  by  ?)  to  be 
the  respective  transport  coefficients.  Of  course  the  calculations  are 
far  more  complicated,  but  bearing  the  above  in  mind  the  calculations 
that  follow  will  seem  clearer. 


We  will  first  put  the  left  hand  side  of  equation  4-39  into  a 
more  convenient  form.  Starting  with  equation  4-38,  we  see  that  in  terms 
of  our  generating  function  we  have 


(4-40) 


n  Ayv  3  o 

U  =  rn  C  - 

o  33. ae  88 

A  y  v 


From  equation  4-40  we  get 


(4-41) 


Ayv  _  Ayv,  r,.Ay  v  .yv  A 
U  =  E,  u  u  u  +  4  (A  u  +  A  u  + 


.vA  y  v 
A  u  ) 


where 


(4-42) 


3,  VB> 

E,  H - 7?  =  4tt  ma  — - - 

Z  p 

c 


K„(3) 

r,  /  3 

~2  =  4tt  ma  - 

c  8 


Thus  the  expression  for  u^V|  can  be  written  as 

o  v 


U  XyV  |  =  (£  +  £  uu.  +  -y  4  +  -j  4  u“i  )u,vu 

o  v  a  2  2  a 


a 


2  a  .  A  y 


,  *y  y  ,  3C  •  \  A  ,  / _  •  A  ,  I A  ,  34  y 

+  (CuM  +  4|M+— yu)u  +(?u  +  4  +  — j  u  )  u 


(4-43) 
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+  (C  +  5  ua,  )AXy  +  C(uXly  +  uUlX) 
I 


Using  the  facts  that. 


a A  a ^  f„p\a  l  °|pn  A|y  y  |  A  ,  u^uV  ,  uXu^ 
A  p  A  a  C.u  +  u  1  )  =  u  1  +  u  1  +  — —  +  — Y~ 

c  c 


and  that 


.Ay  .A  .y  Apa 
A  =  A  A  A 
P  <p 


equation  4-43  becomes 


(4-44)  U  Xuv,  =  WX|J  +  [(?  +  %uw  +  Au  ?I“]UX 
o  v  2  a 

c 


+  [(£  +  \)uX  +  AXa  ?'a]uy  +  [k  +  U  +  )ua,  ]uXuU  + 


.  .  5  a  N  .  Ay 

+  (?  +  3  C  u  |  JA 


a 


where  we’ve  made  use  of  equation  4-17,  and  the  fact' that 


y  Ay  a  r  u 

41  =  AH  c  -  “  o" 

a  2 

c 


y 


Let  us  for  the  moment  consider  the  coefficient  of  the  u  term.  Making 
use  of  equation  3-65  the  coefficient  becomes 


(4-45)  y  _ 
q  = 


a  +  ^§>  (?  + 

- - —  A^  pb  +  AP  C|a=AP  [ - £—p'“+  C,a] 

„  a  a  a  ,  .  P,  ^ 


(y  +  -j) 


(y  3 — 2^ 
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From  (4-42)  we  see  that 


(4-46) 


Ka 

(5  +  =  4tt  ma  (— - p 

c  3  e2 


Using  equation  3-60,  and  equation  4-46  we  see  that 


(4-47) 


(r  +  — ) 

U  2  K  (6) 

c  4  1 


From  the  fact  that  P  = 


(4+^)  V8)  8 

c 

2  K2(8) 

4tt  m  c  a  - —  we  get 

8 


(4-48) 


I*  _ 


pi  =  -  4i  m  c  a 


h. 

B2 


2  K2  |x 

+  4i  m  c  — y  a 1 

6 


and  equation  4-42  we  get. 


(4-49) 


JX  2  ,1  K4.  K3  Jx  2  K3  |X 

4  1  =  4tt  m  c  a  (—  -  — )  —r  6 '  +  4i  m  c  —r  a 1 

3  3  3Z  62 


Substituting  equations  4-47,  48,  49  into  equation  4-45  we  obtain 


(4-50) 


4tt 


2 

m  c  a 


K/  i 
(---)> 
K3  8 


a 


But  we  know  that 


(4-51) 


2 

m  c 

Log  a  =  — —  0  +  const.  (15) 
K 


. 


81  - 


Thus  equation  4-50  becomes 


(4-52) 


2 

m  c 
~K 


U  ~  P 


(5 


(y  i — o) 


-] 


Defining  the  relativistic  enthalpy 


n  to  be, 


we  are  able  to  obtain  the  rather  obvious  identities 


(4-53) 


C  =  n  P 


c 


,  p 

2 

c 


Using  these  identities,  and  the  fact  that 
becomes 


P_ 

2 

c 


P. 

8 


equation  4-52 


(4-54) 


1]  A 


y 


a 


where  we've  used  the  following  results  given  by  Israel  (16); 


(4-55) 


dp  2 

d3  “  n 


and 


Y  =  _  R2  in 

Y-l  d8 


where  y  is  the  ratio 


of  specific  heats,  i.e. 
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Let  us  next  consider  the  coefficient  of 


X  y 
u  u 


in  equation 


4-44,  i.e. 


k  +  a  +  ^f) 

c 


a 


u 


a 


It  can  easily  be  seen  that 


y 


from  which  we  can  write 


(4-56) 


From  equation  4-53  we  can  write 


(4-57) 


C  =  q  P  +  q  P 


Along  the  streamlines  entropy  is  conserved  (17)  i.e.  s  =  0  .  Thus  we 
can  write 


n 


9  n 

8P 

s 


P 


) 


from  which  it  follows  that 


(4-58) 


n 


1 

8 


X 

u 


X 


(18) 
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One  can  also  show  that 


(4-59) 


•  2  2  S 

P  “  P  c  T0  -  9  c  p  | 


=  (P[nB(y-l)  -  y]  -  n  P  B(y-l)}  ua,  .  (19) 


a 


Using  equations  4-57,  58,  59  equation  4-56  becomes 


(4-60) 


5  =  P  +  ^  {-,i^  -  u“|a 

C  1 


Using  equation  4-60,  the  fact  that  (£  i - j)  =  (p  H - )  and  also 

X  C  C 

equation  3-64,  the  coefficient  of  u  u^  becomes, 


Pn  ro  ,  3y,  a  _  a 

— y  [3y  -  5  +  — r]  u  |  =  A0  u  | 

2  q8  a  2  a 

c 


By  a  similar  calculation  the  coefficient  of  the  term 

in  equation  4-44  becomes 


1  r  ,  3y ,  a  _  .  a 

-  —  D  P[3y  -  5  +  — r]  u  |  =  A  u  | 

3  qB  a  3a 


Thus  equation  4-44  can  finally  be  written  as 


, ,  ,,  Xyv  „  Xy  ,  y  X  ,  X  y  .  .  a  XylA  a  .Xy 

(4-61)  U  M  i  =  q  P  a  +  q  u  +  q  u  +  A  u  i  u  u  4-  A  u  i  A 

o  v  z  a  a  a 


where 


84  - 


where 


A  =  -  —  T  3y  -  5  +  -^1 
2  "  2  LJy  n8J 

c 

A  =  -  i  n  p  [3y  -  5  +  ^] 
j  3  n  3 


Next  let  us  consider  the  right  hand  side  of  equation  4-39  in 

^  ^ 

particular  let  us  consider  the  structure  of  X  .  Our  task  now  is 

to  ascertain  the  general  form  of  X^0^  #  From  the  definition  of  Xa^^ 
(see  equation  4-47)  one  can  see  that  there  is  no  special  spatial  direc¬ 
tion  singled  out  .  Thus  we  ask  ourselves  what  is  the  most  general 

form  that  we  can  choose  and  of  what  quantities  will  it  be  composed  of. 
Clearly  the  only  basic  quantities  at  our  disposal  are  u^  ,  and  A^V. 

As  was  stated  earlier  the  quantity  A^V  projects  other  tensors  onto  the 
3-flat  orthogonal  to  u^  .  Thus  the  most  general  combination  of  these 

quantities  which  give  us  a  4th  order  tensor,  and  also  satisfies  the  sym- 
.  j.  va3Xy 

metry  properties  of  X  is 


(4-63) 


xXya8 


A 


X  p  a  8 

u  u  u  u 


+  B  (A 


Xa^y8 


+  A 


ya^Xy 


) 


.  X  pa  8  ,  y.Xa  8 
-I-  G(u  A  u  +  u  A  u 


+ 


X.u8  a 
u  A  u 


+ 


u 


W) 


,  .  Xy  a8  .  w.Xy  a  3  ,  Aa8  X  y. 
+  D  A  A  +  E  (A  u  u  +  A  u  u  ) 


. 
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It  should  be  noted  that  there  are  only  3  independent  components  i.e. 
D  ,  and  E  can  be  expressed  in  terms  of  A  ,  B  ,  and  C  . 

It  can  be  seen  that 


(4-64) 


-8  va8Xy 

A  =  c  u,  u  u  u_  X 
X  y  a  8 


(4-65) 


n  1-4  .  vXyaB 

C  =  —  c  u.  u„A  X 

3  X  8  ya 


To  find  what  B  is  we  first  multiply  equation  4-61  by 


A  A  ,  i.e. 
a8  Xy 


(4-66)  A  A.  Xa3Xy  =  B(AXRA3  +  Ay  A3  )  +  D  Aa3AXyA  A 

a8  Xy  8  X  8  y  a8  Xy 


Equation  4-65  simplifies  to 


(4-67) 


A 

a8  Xy 


xa8Xy 


6B  +  9D  , 


aX  Aa 
since  A  A  , 
a  X 


Similarly  we  can  show  that, 


(4-68) 


A.  A  0  Xa3Xy  =  12B  +  3D  . 
Xa  y8 


Thus  from  equations  4-67,  68  we  see  that. 
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(4-69) 


B  =  ^  (A  A  -  \  A  A  )  x“BX,J 
10  Aa  yg  3  ag  Ay 


Multiplying  equation  4-62  by  u  n  A,  we  get, 

a  g  Ay 


(4-70) 


.  „Ayag  „  4  ^ 

u  u  A,  X  =  3  c  E  . 
a  g  Ay 


But  from  the  definition  of  Xa^y  it  follows  that  g  x^hag  _  q 

Ay 

since  6(P^P^)  =  0  .  So  that  equation  4-70  becomes 


(4-71) 


A  y  agAy  4 

u  u  — x—  X  =  3  c  E 
a  g  2 
c 


But  in  view  of  equation  (4-64)  we  finally  have 


(4-72) 


E  =  i  c2  A 


Similarly  from  equations  4-66,  63  it  follows  that. 


(4-73) 


14  2 

D  =  —  cA--jB  . 


Thus  equations  4-71,  72  show  that  there  are  only  3  independent  components 
in  equation  4-63. 


Our  next  task  is  essentially  to  get  b  in  terms  of  t  , 


ag 


ag 


,a 


and  Q  .  As  was  mentioned  earlier  this  will  be  accomplished  through 


y 


equations  4-22,  23.  Before  doing  this  let  us  first  express  ,  and 


yv 


in  terms  of  b  0  and  moments  of  N  .  Substituting  equation  4-34 
ag  o 
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into  equation  4-5  we  have 


=  ra  c 


N  f  Py  dco  =  m  c  [  N  (C,.  +  C  PV  +  b  0  PaP3)  PU  dm 
Jo  J  o  5  v  ag  * 


which  simplifies  to 


(4-74) 


M/  =  \  b  U  “BlJ  +  i  C  T  WV  +  C5  M  V 
i  ^  ag  o  cvo  5o 

c 


pv 


Similarly,  for  we  get 


(4-75) 


T/V  =  tt  ,  +  1  c  u  aVV  +  uv 

1  4  ag  o  cao  5o 

c 


where  the  fourth  moment  of  N  ,  V  °^V  is  defined  as 

o  o 


(4-76) 


T7  aSyv  _  4 

V  =  m  c 

o 


N  pap3pypv  dm 
o 


But  in  equations  4-74,  75  we  see  that  we  still  have  the  unknowns  , 

and  C<-  to  account  for.  Through  our  fitting  conditions  (i.e.  equations 

4-4,  5)  we  will  be  able  to  express  C..  and  C  u^1  in  terms  of  b  „  . 

5  u  ag 

It  should  be  noted  that  all  we  can  fix  is  u^  .  The  spatial  part  of 

C  remains  arbitrary:  which  is  a  result  of  the  degree  of  arbitrariness 

y 

weTve  left  in  u^  ,  (i.e.  in  not  picking  a  particular  space-time  direc¬ 


tion  for  u  ).  From  one  of  our  fitting  conditions,  equation  4-4,  we  see 


that 


.. 


88  - 


(4-77)  0=M,Pu  =i-b0U“B,Ju  +iC  u  TyV  +  C,  u  M  U  . 

1  y  c^a8o  y  c  v  y  o  5yo 


Substituting  equations  3-58,  62,  and  4-41  into  equation  4-77  we  get. 


0  =  ~  b  u  {5uauSuU  +  3CA(a|Ju6)  +  i  C  u  (yuV  +  PAyV)  +  Cru  (puy) 
2a6y  cvv  5  y 

c 

ct  8 

Assuming  that  g  b  =  0  (20)  ,  and  defining 

Otp 


a  8  _  -v 
ba6U  “  =  H(X)  ’ 


the  above  equation  becomes, 


(4-78) 


C  uV  +  ^  C  = 
v  y  5 


—  (p  +  4  n  p)h(x) 
cy  z 

c 


From  the  other  fitting  condition,  equation  4-5,  we  have 


//  m  yv  1  .  T7  agyv  .  1  _  ayv  , 

(4-79)  0  =  T-  u  u  =  — 0-boV  uu  -I - C  U  uu  + 

v  1  y  v  ^2  a3  o  yvcao  yv 

+  C.  T  UV  u  u 
5  o  y  v 


Using  equations  3-62,  4-41,  and  the  fact  that 


T7  f.  T7  qv 

(4-80)  v“ByV=  4™  +  a“uBuV 

°  8  82  83 

.  2  ,  K4  K3,  ,  (v.y) (a  8) 

+  4tt  m  a  c  { — y  -  — w  4  u  A  u 

3  3 


89  - 


+  4*  m  ac2  ^  {4VVS  +  AVVa  +  a“VV> 

8 

+  4tt  a  m  o2  {^|  -  ^3}  {uauSAPV  +  uVi®6}  , 

e  3 


equation  4-79  becomes 


(4-81) 


K  5K  2K 

0  =  4t  m  c  a  {  — - %  +  — f  }  H(x) 

6  B2  83 

,  3  r  _  a  ,  4  _ 

-fc^Cu-tycC,. 

a  5 


Using  the  recursion  relation 


K 


n+1 


K  ,  +  K 
n-1  8  n 


we  see  that, 


K3  +  !  K4 


K4  K2  +  8  K3  ’ 


thus 


(4-82) 


8K2  48 

S  "  K3  +  ~T  +  ~2  K3 

8 


From  4-42,  3-60,  and  equation  4-82,  equation  4-81  becomes 


(4-83) 


c  ua  +  r  c  c.  = 

a  Ed 


'k  UB  +f)  J2+2Y }  H(x) 

c  c 


(21) 


90  - 


From  equations  4-78  and  4-83  we  are  thus  able  to  solve  for  C  ua  and 

a 


r 


(4-84)  C  ua  = 
a 


qc  (1  Y  )U  +  \)  pc(1_y)  [  (g  +~)t  +  3P] 


+ 


>  H(x)  , 


and 


(4-85) 


C.  = 


u(l-Y)  [(8  +  “)c  +  3P]  ?  °  +  2)(1  Y) 

- _ - § - + - % - 


H  (x) 


where  we  have  also  used  the  fact  that 


(4-86) 


2  r  p 

y  -P  K  = 


c4(1-y) 


.  (22) 


Thus  we  see  that  in  equations  4-74,  75  we  have  ,  and 


yv 


T,  K "  expressed  in  terms  of  b  „  .  Our  next  task  is  to  express  the  Heat 
1  ap 

flow  vector,  and  the  Viscous  S trees  Tensor  in  terms  of  b  n  :  we  do  this 

a(3 

by  means  of  equations  4-22,  23.  Let  us  first  consider  the  heat  flow 
vector,  Qy  .  From  equations  4-74,  3-58,  3-62,  4-41,  the  first  term  on 
the  right  hand  side  of  equation  4-22  becomes 


(4-87)  -  nc2  Miy  =  -  q  £H(x)uP  -  2q  ?  b  A^u3  -  ^  H(x)uy 

c 


-  q  yc(C  uV)uy  -  P  q  c  C  AyV  -  q  p  c^  Ch  uy 
V  v  5 


Upon  substituting  equations  4-84,  85  into  equation  4-87  we  get 


-  91 


( I- Y )  (5  +  ^y) 

-  n  c2  M1y  =  {-  n  ^  +  n  c4 - - - —  (y2  -  pO>  h(x)  uy 

c  P 

-  2  t,  b  Aay  u3  -  P  q  c  C  AyV 
a3  v 

By  means  of  equation  4-86,  the  above  simplifies  to 

(4-88)  -  n  c"  M  11  =  -  2  ri?  b  Aay  -  P  n  c  C  AyV 

1  ap  v 


Similarly,  by  substituting  equations  3-62,  4-41,  80  into  equation  4-74, 
the  second  terms  of  equation  4-22  becomes 


(4-89)  -  Ay  T  y  °  u  =  4tt  m  a  c2(— y  - 

P  1  V  b2 

Thus  by  adding  equations  4-87,  and  4-89 


^)(2)  b«8  a1J“uB  +  s  c  c« 
p 

,  equation  4-22  becomes 


A 


ya 


(4-90)  0U  =  -  2  P  [q2  -  ^  -  l]b  Aay  u3 

3  a3 

and  using  the  fact  that 


Y 


Y- 


1  (23) 


equation  4-97  finally  yields 


2P 

B2 


( 


y^I)ba3 


ay 


3 

u 


(4-91) 


92  - 


In  order  to  replace  the  quantity 
thermodynamic  variables  we  have  to  go  back 
use  of  equations  4-61,  62,  63,  71,  i.e. 


b 

to 


ay  6 
aB  4  U 

equation 


by  gradients  of 
4-39,  and  make 


n,  qo\  -n  Xy  ,  y  A  Ay,.  a  A  y  ,  A  a  .Ay 

(.4-92)  qPa  +  q  u  3-qu  +  A  u  i  u  u  +  u  i  A 


2  a 


3  a 


,  r »  A  y  a  3  ,  ,,,AAaAy8  ,  .ya.A3N  ,  A  ya  3  , 
=  -b  {A  u  u  u  u  +  B (A  A  +  A  A  )  +  C(u  u  u  + 


.  y.Aa  3  ,  A.y8  a  ,  y.Ay  aN  ,  „  .Ay  .a8 
+  u  A  u  +  u  A  u  +  u  A  u)+DA  A 


+ 


,  w.Ay  a  8  ,  Aa8  A  y.  , 
+  E(A  u  u  -I-  A  u  u  )} 


Multiplying  equation  4-92  by 


u.  u  we  get , 
Ay  6  5 


(4-93) 


4  . 

c  A2  u 


a 


a 


=  -  c 


A  b 


a3 


a  3 

u  u 


-  E 


2  K 

c  b 


a8 


a  3 

u  u 


where  we  have  used  the  fact  that  ua^  =uq^=0,  (this  follows  easily 

y  y 

from  definition  of  ,  and  q^) .  But  in  view  of  equation  4-72,  equa¬ 

tion  4-93  yields  the  result. 


(4-94) 


Multiplying  equation 
we  get 


4  A 

3  a 

u  I 

a 


H  (x)  . 


4-92  by  A 

Ay 


and  noting  the  fact  that 


93  - 


(4-95) 


O  A  a 

3  A_  u  | 

3  a 


-  0^|  +  ^~|)  H(x)  -  3E  H(x) 
c  c 


But  in  view  of  equation  4-73,  equation  4-95  simplifies  to, 


(4-96) 


a  4  A  c' 

A3  =  -  9  — 
u 


a 


Multiplying  equation  4-92  by  u  ,  and  noting  equation  4-94,  we  get 

r* 


(4—97 ) 


4  o  r1  ru  .  X  Ot  3  \ 

q  =  -  2C  (b  ^  A  u) 
aS 


Thus  substituting  4-97  into  equation  4-91,  and  noting  the  definition 


of  q^  (i.e.  equation  4-62)  we  finally  get 


(4-98) 


*4  =  h^)]2  AV“ 

p 


thus  the  thermal  conductivity  X  is 


(4-99) 


X  =  F—  (—■  )  1  ^ 

nc  L32  (y-i 


We  still  have  to  calculate  C  as  defined  in  equation  4-65.  This  will 
be  done  at  a  later  stage. 


Let  us  next  consider  the  viscous  stress  tensor.  From  equations 


4-34,  75,  80,  41,  and  3-62  we  have 


94  - 


(4-100) 


t  „  =  ~  b  { 4t t  mac^  — |  (2AY  A^  +  A^A  ) 

^3  a  3  a3 


aS  2  y6 
c 


?  KA  A  1 

+  4tt  mac  (— r- - qO  uYu  A  +  —  r,  C  uY)  A  0 

^2  ^3  a3  c  y  a8 


+  Cc  P  A  .  . 

5  a3 


In  view  of  equations  3-60,  4-42,  82,  equation  4-100  becomes 


(4-101) 


Lb  ■  f  (2  VW  +  (¥  H  +  V  (p  +  ¥  H  + 

c  6  c 


+  -  C  C  uY  +  C..P)  A 
c  y  5  a6 


6  6 

In  order  to  eliminate  the  quantity  b  „A  A  „  we  once  again 

Y  <S  a  6 

go  back  to  equation  4-92.  Substituting  equations  4-94,  95,  96,  103,  104 
into  equation  4-91,  and  noting  equations  4-72,  73,  81  we  have 


(4-102)  q  P  a 


Xy  -  \  \  H(x) 
c 


-  2  b  b  y  vB 

a8 


Making  use  of  equation  4-102,  equation  4-101  becomes 


(4-103)  t 


a8  8 


(dP)  _  f5  q  P  1  .  5qP . 

- r  °aB  +  [3  TV  H  +  ~2  (P  +  “T5  H 

6c  c 


+  -  c  (C  uY)  +  CP]  A  . 
c  y  5  a3 


Let  us  now  calculate  the  coefficient  of  A  :  Upon  substitu- 

aP 

ting  in  the  values  of  C^uY  an<3  C5  ^rom  equations  4-84,  85,  93  we  get 


95  - 


{T  t  ~2  +  "T  +  -  PC)(C  +  -%)  +  (PC  -  yP)  [  (3  +  ~  K  +  3P  ]  ] 

c  c  c  B 


(1-y) 


P2  p2 

Using  the  identities  (24)  pz  -  nP  =  —  ,  and  pc  -  JP  =  -h- 

c  3c 

the  above  simplifies  to 


[UB  -  £.] 


{ 


20  n  P 

3  Q  2 

3c 


4P  (  2  5n  ..  4yP  ,  2  5q 

T  (n  -  T  '  x)  +  2  (n  -  T 


-  i) 


4qyP 

3c2 


}  H  (x) 


Making  use  of  equation  4-55  the  above  yields 


(4-104) 


-'V  [3y  -  5  +  -fq  H(x) 

3cZB  Bn 


But  by  equation  4-94,  4-104  becomes 


(4-105) 


1  ,Pru2  1  n2  a 

-  A  (T>  u  let 

c  1 


3y 

where  0,  =  (3y  -  5  +  .  Expression  4-105  is  the  coefficient  of 

3n 

A  .  in  equation  4-103.  Therefore  the  viscous  stress  tensor  becomes, 
a3 


(4-106) 


a3 


r(nP)  ,  ri  ,Ptk2  1  _2,  X  . 

[  3B  ^  a3  A  2  3  -1  U  |x  Aa3 

c  1 


But  this  is  just  the  phenomenological  representation  of  the 
viscous  stress  tensor  (see  equation  4-17).  Thus  we  are  able  to  identify 
the  shear  viscosity  and  bulk  viscosity  with  the  coefficients  of  , 


96 


and  u  | ,  A  „ 
|  A  a3 


respectively,  i.e. 


(4-107) 


v  = 


(TIP)' 

3B 


,Pr)fh 2 

K  =  (— — ) 


A3 


As  was  mentioned  in  the  case  of  conductivity  the  quantities 
A  ,  B  ,  C  have  to  be  calculated  before  we  have  the  coefficients  just 
as  functions  of  temperature  T  (or  3) ,  and  y  .  This  next  section  will 
involve  choosing  proper  coordinates,  basis  vectors,  and  then  integrating 
to  get  A  ,  B  ,  C  . 


Integrations 


Preliminaries ♦ 

In  order  to  perform  the  integrations  we  choose  proper  coordi¬ 
nates  with  the  time  axis  along  u  ,  so  that  u  =  (0,0,0,-c)  .  In  order 

y  U 

to  specify  the  relative  momentum  of  the  incoming  particles  we  choose  an 

orthonormal  triad  e,  ,  e0  ,  and  e„  in  the  centroidal  3-flat.  The 

~1  ~z  ~3 

comoonents  of  the  vectors  e.  ,  e0  ,  e0  are,  (see  figure  1) 
r  ~1  ~2  ~3  . 


(4-108)  &C/o\  ~  (cosh  x  sin  6  cos  <j>>  cosh  x  sin  6  sic  $  >  cosh  x  cos  0, 


(3) 


sinh  x)  > 


ct 

e^  =  (sin  <J>,  cos  <f> ,  0,0)  , 


a 


(1) 


=  (sin  0  cos  <J) ,  sin  0  sin  <f> ,  cos  0,  0) 
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It  can  easily  be  seen  that  do  indeed  form  an  orthonormal 

triad,  i.e. 


(4-109) 


~(a)  ~(b)  Yag  G(a)  e(b)  ^ab  ’ 


where  y  is  the  Minkowski  metric,  i.e.  diag.  (1,1, 1,-1)  ,  and 
a,3>  =  1,2, 3, 4  ,  and  a,b,  =  1,2,3  . 


It  can  also  be  shown  that  e,  ,  en ,  e„  due  in  fact  lie  in  the 

~1  ~2  ~3 


centroidal  3-flat; 


e,  .  •P=ye/NP 
~(a)  ~  a6  (a) 


“  =  0  V 


e.q.  consider  the  case  for  a  =  3  ,  we  know  that 


^=-a 


P  =  m  c(sinh  x  sin  9  cos  sinh  x  sin  9  sin  4>  >  sinh  x 

cos  0,  cosh  x] 


thus 


—  a  — 8  —  —  2  2 

e„  •  P  =  y  n  e/nN  P  =  m  c(sinh  x  cosh  x  sin  0  cos  d>  + 
~3  ~  a3  (3) 


—  —  2  2 
+  sinh  x  cosh  x  sin  9  sin  <j>  + 


+  sinh  x  cosh  x  cos  "0  -  cosh  x 


x  sinh  x)  =  m  c  (sinh  x  cosh  x 


x 
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FIGURE  1 


99  - 


x  [sin  0  (cos  <j)  +  sin  <f>)  +  cos  0]  -  cosh  x  sinh 
=  0  . 


Thus  using  polar  coordinates  'O  t  and  e  we  can  specify  the 


•  ct 

normalized  relative  momentum  vector  g  /g  by 


a 


(4-110)  =  sin^  cos  e  e^  N  4-  sln^O  sin  e  +  cos^  cos°Lx 

g  (1)  (2)  (3) 


where  *0  and  e  are  measured  as  indicated  by  figure  2.  Using  equation 

ot 

4-110,  it  can  be  shown  that  g  /g  as  defined  by  equation  4-110  does 
indeed  satisfy  the  condition. 


a  3 

Y  =  1 

g  g  ct3  g  g 


a  g 

.  fa  = 


In  order  to  specify  the  normalized  relative  momentum  of  the 
outgoing  particles,  i.e.  g*/g  ,  we  make  use  of  another  orthonormal 
triad  E,  ,  E„  ,  and  E_  ,  (see  figure  3) ,  in  the  centroidal  3-flat. 
The  vectors  in  the  triad  are  defined  as 


a 


Ea  =  &- 
(3)  g  ’ 


,a  a  ,  a 

"(2)  =  "  Sln  £  e(1)  +  cos  €  e(2)  ’ 


„a  S)  a  ,  S)  a  .  f)  a 

E^)  =  cos  'V  cos  e  eQ)  +  cos'V  Sm  e  e^)  -  sm'V  e(3)  » 


where  V  ,€ 


a , 

are  those  values  that  specify  g  /g  . 


' 
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FIGURE  2 


FIGURE  3 


FIGURE  4 
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Thus  relative  to  the  triad  E1  ,  E~  ,  and  E_  ,  the  normal- 
ized  relative  momentum  of  the  outgoing  particles  is  given  by 

§  & 

(4-112)  —  =  sin  0  cos  $  E.,  +  sin  0  sin  $  E„  +  cos  0  E„ 

g  ~1  ~2  ~3 

where  0  ,  and  $  are  measured  as  indicated  by  figure  4. 

Now  that  some  of  the  preliminaries  have  been  discussed  let  us 
now  calculate  the  coefficients  A  ,  B  and  C  .  Let  us  first  consider  A  . 
From  equation  4-64  and  the  definition  of  we  see  that 


A  = 


1  -5 

—  me 
4 


/*  >» 


J  J  J  J 


W  N 

o 


’N  u  u06{PaP^}  u  u  dfP^P^ldm  d'w  da)  d'a) 
o  a  6  X  p 


9 


which  becomes 


(4-113) 


A  1  ~5 

A  =  7  me 
4 


W  N  ?N 


o 


{6[E2]}2 


da)  d '  a)  dm  d 1 


* 


0) 


where 


(4-114) 


«[EZ]  =  ku^p11)2  +  (U|j'pm)z  -  (uup"m)z  - 


,U,2 


* 


But  we  know  from  equations  2-25,  27  that 


P 

U 


2 
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U 


n  ,  1  2,2  -*  1  * 

(1  +  4  8  >  P  y  -  2  8  y 


(4-115) 


fP 


'P 


n  ,  1  2,2  —  1  * 

y  =  (1  +  4  g  >  Py  +  2  8  y 

1 

*  1  2,2—*  1  * 

=  (1  +  4  8  >  p  U  +  2  8  y 


Substituting  equation  4-11  into  equation  4-113,  and  using  the  fact 

—  — k 

that  P  =  P  ,  we  get 

y  y 


6[*2]  =  |  t(u SV)2  -  (g*Uu  )2]  . 

^  M  M 


But  since  we  have  u^  =  (0,0,0,-c)  ,  we  finally  get 


(4-116) 


2  1  2  4  2  *4  2 

\  c  [(gV  -  (g  V] 


But  inview  of  the  preliminaries  just  discussed  we  know  that 


(4-117) 


4 
(3) 


g  =  g  e,^,  cos-C?  =  g  cosyO  sinh  x  » 


*4  0  4  4  0 

g  =  -g  sin'V  sin  ©  cos  $  +  g  e(3)  cos  </  cos  0 

(4-118)  g  ^  =  g  sinh  x  (cosd9  cos  0  -  sin^  sin  0  cos  $) 


Thus  from  equations  4-116,  117,  118,  125,  126  we  can  now  write 


22  144  4  —  r  2  n  4  4o  4  4 

(4-119)  [&{&}]  =  c  (g)  sinh  x  (cos  *\)  sin  0  +  sin  a/  sin  0  cos  $ 
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+  4  sin  2d  cos  2d  sin^O  cos2 


2  2  O  2  O  A  2 

in“0  cos~0  cos  $  -  2  sin  a/  cos  a/  sin  0  cos  $ 


-  4  sin2,v)  cos sin^G  cos  0  cos^$  +  4  sinh^  cos  sin^G  cos  Q  cos  $  } 


Before  performing  the  integrations,  let  us  consider  the  other 

terms  in  equation  4-113.  The  term  N  fN  becomes 

o  o 


N  ’N 
o  o 


2  —  (Py  +  TPU)u 

a  e  c 


y  = 


2  r28  n  ,  1  2. 2  —  y  , 

a  exp  [—  (1  +  g  )  u^  P  ] 


which  simplifies  to 


(4-120) 


N  'N 
o  o 


1  2  1/2  — 
2  "26(1  +  g  )  7  cosh  x 


a  e 


Making  use  of  equations  2-37,  38,  we  have 


do)  d’w  =  (1  +  g2)2  g2  doj  dg  sin  V  d  de 


ini?  d^ 


(4-121)  do)  d'm 


=  (1  4-  g2)2  g2  sin^  sinh^  x  sin  9  d  g  d  x  dGd^d^de 


and 


(4-122)  dm*  dfo>* 


(1  +^(gV)2  (gV  cto*  dg*  sin  0  dG  d$ 


From  equation  2-35  we  know  that 


W(P,  rP:P*’P*)  =  (1  +  ^  g2)  2g  1  a( g,0)6(jj(P  ^ 


P^)6(g  -g) 
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Thus  making  use  of  equations  4-120,  121,  122  and  equation  2-35,  integra 

— : k  ■& 

ting  over  dco  dg  ,  and  making  use  of  the  delta  functions,  equation 
4-113  becomes 


(4-123) 


a  1  2-5 

A  =  —  m  a  c 
4 


.2,-2  3 


a(g,0)  exp  [-23(1  +  g2)2  cosh  x]  x 


1 

1  2,2 


2  - 


x  [<${£'*"}]*’  g"’  (1  +^gz')i’  sinh^  x  sin  6  sin^  sin  0 


x  dg  dx  d0  d<j>  d“$de  d0  d$ 

Substituting  equation  4-121  into  equation  4-123,  and  integrating  over 
0,cJ>,  , £  ,^and  $  we  obtain 


a  16  3  2  -1 

A  =  -qr-  tt  a  m  c 


1 

sinh^  x  exP  [“23(1  +  ^  g2)2  cosh  x] 


x  (—  sin^0  +  sin20  cos20)  (1  +  ~  g2)2  g  a(g,0)dg  d0  dx 


1  2  1/2 

Letting  y  =  23(1  +  t  g  )  the  above  becomes 


(4-124) 


*  16  32  -1 

A  =  — T-  TT  a  C 


r  r 


fy  C°Sh  X  sinh6  x(^)  g7  * 


x  ad  sink  +  sink  cosh)  dg  d©  dx 


Using  equation  3-55,  equation  4-124  becomes  on  integrating  over  x 


(4-125) 


A  = 


3  2-1 

8tt  ma  c 


,  K  (y) 

sin  0  — = —  g  o(g,0)  dg  d0 
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Let  us  next  consider  the  coefficient  B  (see  equation  4-69). 
There  are  two  integrations  that  have  to  be  performed.  Breaking  up  the 
calculation  of  B  into  two  different  integrations  we  have. 


(4-126) 


ot3 


)  xaBXlJ 


We  have  essentially  already  calculated  G  in  our  calculation 
of  A  .  Using  identity 


A 

yv 


pypv 


(u  Py)2  „2 

_ y _ i  =  — 

2  2 
c  c 


1 


we  see  that 


A«b\u  X 


a(3Ay 


1 

4  m 


W  N  '  N  [  6  }  ]  2 

o  o 


dco  d 1  oa  dm  d  ’ 


0) 


since 


(u  PM)2  (u  ?V)2  (u  Pu) 

«[— i C - 1]  =  «[— h 1  -  «[1]  -  «[— 1 j— 

c  c  c 


2 

-3 


But  we’ve  already  calculated  the  preceding  integral,  so  that  it  can  easily 
be  seen  that 


(4-127) 


A  .A,  X 
a$  Xy 


a3Xp 


Q  3  23 

ott  ma  c 


o  Ko(y) 
sin  0  — x — 


g  a(g,0)  dg  d0 


3 
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In  calculating  F  we  have  to  calculate 


(4-128)  A  A.  Xa6Xy  =  y  m  c3 
aX  8ij  4 


W  N  'N  A  ,Ad  6{paP3}  6{PXPy} 
o  o  aX  6y 


k  k 

dw  d’w  dco  d'w 


But  since  in  proper  coordinates  A^  ,  and  A  have  components  (1,1, 1,0) 


then  4-128  becomes 


(4-129)  A  Aq  Xa3Xy  =  f  m  c3  [• 
aX  8y  4  I 


W  N  ’N  6{P  P,}  6{Papb} 
o  o  a  b 


k  k 

doo  d’to  doo  d’m 


where  a,b  =  1,2,3  . 


— 3.  — b  —  — 

Using  equation  4-115  and  the  fact  that  P  =  P  ,  and  P^  =  P^ 


we  can  write 


(4-130) 


6 [PaPb ]  [gagb 

6[paPb^  =  2  [gagb 


*a  *b, 
g  g  J 


.8  J 


Thus  we  can  write 


(4-131)  6 [ P ap b 1 6 CPaPb]  =  \  ((ga g3)  +  (8  ag  a)  “  2 (§aS  a)2^  • 


Since 


a 

g  g 


a 


=  g 


*  a  * 


a 


=  g 


we  can  write 


- 
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(4-132) 


a  r  2  r  4v2-> 

g  ga  =  [g  +  (g  )  } 


*a  *  r  2  ,  *4.2, 

g  g  „  =  (g  +  (g  )  > 


It  can  also  be  easily  seen  that 


>a 


a 


=  g  g 


ra 


(3)a 


=  g  cos  0 


Thus  we  can  write 


//  *a  2  *44 

(4-133)  g  g  =  g  cos  0  +  g  g 

3 

Making  use  of  coordinate  system  we  find  that 


g  =  g  sinh  x  cos 


-0 


(4-134) 


g  ^  =  -g  sinh  x  sin^  sin  0  cos  $  +  g  sinh  x  cos^  cos  0 
g3g  =  g^(l  +  sinh^  x  cos^) 

3 

*a  *  2  ,  .  2  -  r  .  2  Q  .  2n  2.  ,  2  Q  2n 

g  g  =  g  11  +  sinh  x  [sm  sm  0  cos  $  4-  cos  V  cos  0 

3 

-  2  sin^  cos  $  sin  0  cos  0  cos  $]} 


g  ag  =  g^  {cos  0  +  sinh^  x  [cos  2J  cos  0  -  sin'$  cos  ^  sin  0  cos$>]} 


Setting  up  our  integrations  in  exactly  the  same  fashion  as  we  did  in  calcula¬ 
ting,  and  after  many  integrations  over  X»9>(j)>  ,£,^and  $  we  finally  arrive 


at 
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(4-135)  A  A  X 
aA  8p 


oc8Xp 


n  3  2  3 

oft  ma  c 


^2  ^^3  3  7 

(Kx(y)  +  — -  +  — 2)  sin  0  g 

y  y 


o(g,0)  dg  d0  , 


—  12  1/2 

where  in  integrating  over  x  >  we  once  again  let  y  =  28(1  +  g  ) 

and  made  use  of  equation  3-55,  Thus  from  equations  4-120,  121,  129  the 

coefficient  B  (see  equation  4-69)  becomes 


,323 

(4-136)  S4nip. 

0  8 


2K  (y)  K  (y) 

(K^y)  +  — -  +  — 2 — )  sin  0  § 


o(g , 0)  dg  d0  o 


From  equations  4-65,  39,  the  coefficient  C  is 


(4-137)  C  =  m  c  1 


W  N  'N  u  A  u  6{PaPB}6{PAPy}dw  d'co  dm*  d* 
o  o  X  pa  8 


w 


In  proper  coordinates  equation  4-137  reduces  to 


W  N  'N  [6{PT  }]  [6{P4Pa>]  dw  d'co  dw  d'w  , 

00  a 


where  a  =  1,2,3. 

From  equation  4-130  it  follows  that 


(4-138)  C  =  m  c 


a  14  *4  * 

6{P  Pa>  =  2  (g  8a  "  g  g  a} 


.  r_4_,a,  1,4a  *4  *a. 

6{P  P  }  =  1  (g  g  -  8  g  ) 


(4-139) 
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Making  use  of  equation  4-139,  134  and  once  again  embarking  upon  many 
integrations  we  finally  arrive  at. 


(4-140) 


C 


2  3  2 

—  name 


2K^  2  7 

( -  +  —x)  sin  0  g  a  (g  ,0)  dg  d0  . 

y  y 
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CHAPTER  5 


This  concluding  chapter  will  be  devoted  to  a  brief  discussion 
of  the  results,  the  limitations,  and  the  possible  extensions  of  the 
theory  of  the  previous  chapters. 

If  we  compute  the  asymptotic  behavior  of  the  bulk  viscosity, 
using  the  necessary  expressions  given  in  chapter  4  we  find  that  for  the 
classical  (3-*30)  and  the  ultrarelativistic  (p-K))  cases  the  bulk  vis¬ 
cosity  vanishes  for  a  simple  gas  (1).  We  also  see  why  it  is  impossible 
to  have  an  isentropic  expansion  of  a  relativistic  gas  (2). 

As  an  application  of  the  above  let  us  consider  the  problem  of 
a  spherically  symmetric  expanding  universe  filled  with  matter,  i.e.  the 
Robertson-Walker  Universe.  In  the  usual  discussions  that  were  given 
about  the  Robertson-Walker  Universe,  we  had  a  universe  filled  with  some 
fluid  that  had  no  shearing  or  heat  conduction,  and  that  is  expanding. 

This  would  imply  that  it  is  expanding  with  a  constant  entropy  (i.e.  no 
dissipative  agents  are  present).  But  from  our  previous  considerations 
(see  equation  3-49)  we  know  that  the  equilibrium  state  given  by  a  constant 
entropy,  i.e.  |  =  0  ,  implies  the  existence  of  a  time-like  Killing 
vector  field.  But  it  is  known  that  a  Robertson-Walker  type  Universe  does 
not  admit  a  time-like  Killing  vector  field.  Thus  it  seems  as  if  we  have 
some  sort  of  a  paradox.  But  in  actuality  there  is  no  paradox  because  of 
the  fact  that  even  though  we  don’t  have  any  shearing  or  heat  conduction  we 
still  have  bulk  viscosity,  and  thus  we  have  a  dissipative  agent.  To  see 
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this  more  explicitly,  let  us  consider  our  universe  in  question  to  be 
described  by  the  Robertson-Walker  metric, 

(5-1)  ds2  =  dt2  -  R2 (t)  +  r2  d02  +  r2  sin20  d<j)2]  . 

1-fcr 

Since  the  above  line  element  is  written  in  the  comoving  coordinate 

12  3  4 

system,  we  thus  have  u  =  u  =  u  =  0  ,  and  u  =1  for  velocity  u  . 
Since  we  have  no  shearing  or  heat  conduction,  we  see  that  in  the  comoving 
frame  the  Energy-Momentum  tensor  given  in  equation  4-13  becomes 


,yv  _  yv  yv  X 

=  T  -  k  gH  u 
o 


y  v  ,  X 

=  y  u  u  +  (P-ku 


)  A 


yv 


From  which  we  deduce  that 


(5-2) 


and 


(5-3)  Tlj  =  (P-kuA,  )g1J  =  (P-<-r“  )g1J 

A  *4CX 


where  i,j  =  1,2,3  . 

In  view  of  equation  5-1,  equation  5-3  yields 


(5-4) 


T1J  =  (P  -  3k|)  glj 


But  we  know  that  in  the  early  "hot"  stage  of  the  universe, 
and  the  later  "cold"  stage  of  the  universe  k  vanishes.  So  that  in 


-  112 


yv 

these  stages  T  ,  as  expressed  by  equation  5-4,  is  just  the  Energy- 
Momentum  tensor  of  a  perfect  fluid.  But  in  the  intermediary  stages  k 
is  non-vanishing  and  hence  T^V  is  an  imperfect  fluid,  which  in  turn 
means  entropy  production. 

As  a  result  of  the  above  we  are  led  to  the  question  of  Helium 

production.  We  know  that  Helium  production  takes  place  at  temperatures 

9 

of  the  order  of  10  K  .  This  temperature  is  intermediate  between  the 
ultrarelativistic  and  the  classical  ranges.  But  if  one  examines  equation 
5-4  it  can  be  seen  that  the  bulk  viscosity  opposes  the  pressure.  The 
question  then  to  ask  is:  MCan  the  bulk  viscosity  counter  act  the  pres¬ 
sure  enough  to  slow  down  the  expansion  rate  so  that  it  spends  a  slightly 

9 

longer  time  around  T  ~  10  K  and  thus  produces  a  significant  amount  of 
more  helium  than  would  have  been  produced  otherwise?".  There  are  essen¬ 
tially  two  main  reasons  why  one  would  expect  the  answer  to  be  no.  Firstly 
the  pressures  are  so  large,  and  the  expansion  rate  (R/ R)  so  small  that 
bulk  viscosity  couldn’t  affect  the  expansion  to  any  appreciable  extent. 

The  second  reason  is  that  the  measured  abundance  of  Helium  is  in  close 
agreement  with  theoretical  predictions  which  neglect  bulk  viscosity  (3). 
One  should  note  that  this  second  reason  given  is  not  so  devastating.  The 
reason  being  that  these  measured  abundances  are  gotten  from  "new"  star 
populations.  Accurate  measurements  have  not  yet  been  made  for  the  old 
stars  in  the  galactic  halo.  The  abundance  in  these  old  stars  is  closer 
to  the  amount  of  helium  produced  cosmologically . 

Another  interesting  aspect  of  bulk  viscosity  is  whether  it  can 
account  for  the  present  entropy  per  baryon.  Weinberg  (4)  has  shown  by 
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rather  general  arguments  that  bulk  viscosity  could  not  have  given  rise 
to  the  present  entropy  per  Baryon  for  the  present  microwave  background, 
in  a  Robertson-Walker  type  universe.  What  he  essentially  does  is  to 
first  compute  the  rate  of  increase  of  the  entropy  per  particle  through 
the  thermodynamical  relation 


(5-5) 


T  ds  =  du  +  P  d  v  , 


where  v  is  the  comoving  volume,  s  ,  and  u  are  the  total  entropy, 
and  the  total  energy  respectively  in  the  comoving  frame.  From  equation 
5-5  we  get 


(5-6) 


T  s  =  u  +  P  v 


where  the  dot  denotes  the  derivative  with  respect  to  time.  By  means  of 
Einsteins  field  equations  one  can  prove  that 


(5-7) 


~  (yR3)  =  -3  R2R  (P  -  3k  |) 


and 


(5-8) 


(nR3)  =  0  , 


where  n  is  the  density  of  particles. 


If  a  is  the  entropy  per  particle  then  the  total  energy  s 


is 
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(5-9) 


s 


n  k  a  R 


The  total  energy  u  is 


(5-10)  u  =  y  R3 

and  the  volume  v  is 


(5-11) 


In  view  of  equations  5-7,  8,  9,  10,  11  equation  5-6  becomes 


(5-12) 


a 


•  2 

9  <  R 
n  k  T  R2 


If  the  bulk  viscosity  arises  from  some  form  of  radiation  then  to  first 
order  in  x  ,  the  mean  free  time,  equation  5-12  reduces  to  (5), 


(5-13)  i  -  If  1  . .  .  . 

adiabatic 

where  P  is  the  adiabatic  rate  of  change  of  the  radiation  entropy;  i.e. 


(5-14) 


P  =  Hr  ( 


4bT' 


dt  3  nx 


)] 


adiabatic 


Then  by  some  general  arguments  he  establishes  limits  to 
The  result  he  gets  is  that  (6) 


a 


(5-15) 


o  i  I  p|  . 
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i.e.  "the  rate  of  entropy  will  always  be  less  than  if  the  radiation  or 
matter  gained  entropy  at  the  rate  |P|  and  the  matter  or  radiation  had 
fixed  entropy".  (7).  Weinberg  thus  arrives  at  two  conditions  that  must 
be  satisfied  if  there  is  to  be  an  appreaciable  production  of  entropy  in 
a  homogeneous  isotropic  expansion: 


(a)  "The  matter  and  radiation  must  both  contribute  an  appreci¬ 
able  share  of  the  total  energy.  If  either  the  radiation 
entropy  or  the  matter  entropy  contributes  only  a  small 
fraction  of  the  total  entropy  throughout  some  time  interval, 
then  in  an  adiabatic  expansion  the  change  of  radiation 
entropy  would  have  to  be  a  small  fraction  of  the  total 
entropy,  so  that  the  growth  of  entropy  could  at  most  be  a 
small  fraction  of  the  total  entropy"  (8). 

(b)  We  see  that  in  view  of  equation  5-15  equation  5-13  implies 
that  for  maximum  entropy  growth,  we  must  have  (9), 


Weinberg  then  shows  that  for  various  kinds  of  radiations  i.e. 
photons,  electron-type  neutrinos,  muon-neutrinos,  and  gravitons,  the  above 
conditions  aren’t  met.  Which  in  turn  leads  Weinberg  to  the  result  that 
bulk  viscosity  could  not  account  for  the  large  entropy  per  baryon  of  the 
present  microwave  background. 
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For  a  detailed  verification  of  Weinbergs  conclusions,  one 
would  have  to  approach  the  problem  from  a  kinetic  theory  point  of  view. 
Clearly  the  kinetic  theory  formalism  as  it  is  discussed  in  this  thesis 
would  be  to  naive  to  describe  the  above  problem.  One  would  have  to  formu¬ 
late  a  theory  for  mixtures,  i.e.  allowing  different  species  of  particles. 

In  formulating  this  theory  for  mixtures  one  would  to  allow  for  more  general 
forms  of  interactions  than  close  binary  collisions.  One  would  have  to 
allow  collisions  in  which  particle  number  is  not  conserved.  One  would  have 
to  also  allow  multiple  collisions  (more  than  binary).  There  are  also  two 
processes  that  must  be  accounted  for  if  the  theory  is  to  be  a  physically 
meaningful  one:  the  first  is  Bremsstrahlung ,  and  the  second  is  the  crea¬ 
tion  and  annihilation  of  particles. 

Finally  we  come  to  the  question  of  the  validity  of  the  Relativis¬ 
tic  Boltzmann  equation.  In  the  case  of  the  non- relativistic  Boltzmann 
equation,  it  can  be  shown  that  it  is  derivable  from  the  more  basic  Liouville 
equation  after  a  series  of  approximations.  Some  authors  (10)  have  questioned 
the  validity  of  the  Relativistic  Boltzmann  equation  because  it  has  yet  to 
be  derived  from  the  more  basic  equations.  But  this  in  turn  raises  the  very 
fascinating  question  as  to  whether  one  can  formulate  a  Relativistic  version 
of  Liouvilles  theorem  for  interacting  particles,  and  thus  give  rise  to  a 


Relativistic  Statistical  Mechanics. 


. 
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FOOTNOTES 
Chapter  I 


(1)  Of  course  we  must  employ  the  usual  assumption  about  d  rd  v  , 

i.e.  that  it  is  large  enough  to  contain  a  large  number  of  particles, 
and  yet  small  enough  to  be  insignificant  with  respect  to  macroscopic 
dimensions.  We  also  assume  that  the  number  of  particles  is  large 
enough,  and  the  variation  of  N  is  not  great  between  successive 
elemental  volumes,  so  that  we  can  replace  everything  by  continuous 
functions,  and  summations  by  integrations. 

(2)  Actually  we  are  looking  for  quite  a  chaotic  or  unstructured  system. 
In  fact  the  most  chaoticl 

(3)  We  are  of  course  assuming  that  the  cells  all  have  an  equal  apriori 
probability  of  being  filled  by  any  one  particle, 

(4)  See  Tolman  "Statistical  Mechanics". 

(5)  See  J.L.  Synge,  "Relativity:  The  Special  Theory",  North  Holland 
Publishing  Company,  1965,  p.  67. 


(6) 


See  Appendix  A. 
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(7) 


i  r  i 

t  n 

=  s  1  t 

n  | 

• 

1  r  1 

1  r 

r 1 

"This  formula  connects  the  3- volumes  s,s  of  cross  sections 

r  — r 

normal  to  n  ,  n  respectively,  of  a  world  tube,  t  being  any 
tangent  vector  to  the  tube." 

See  J.L.  Synge  "Relativity:  The  Special  Theory",  Chapter  8,  section 
6,  North  Holland  Publishing  Company,  1965. 


(8)  W.  Israel,  J.  of  Math.  Phys .  1963  #9,  p.  1164, 

(9)  J.L.  Synge  "The  Relativistic  Gas",  Chapter  4,  North  Holland  Publish¬ 

ing  Company. 

Chapter  II 

(1)  Liouvilles  Theorem. 

(2)  See  p.  2-2. 

(3)  For  proof  see  W.  Israel,  J,  of  Math.  Phys.,  Vol,  4,  #9,  p.  1180. 

(4)  For  proof  see  W.  Israel,  J.  of  Math.  Phys.,  Vol.  4,  #9,  p.  1180. 

(5)  It  should  be  noted  that  the  expression  for  the  scattering  cross 
section,  a(g,9)  ,  that  appears  in  equation  2-41  differs  from  the 
conventional  scattering  cross  section  given  in  physics.  The  reason 
for  this  lies  in  what  we  choose  to  call  the  relative  3-velocity  in 
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dxi 

the  center  of  mass  frame.  The  usual  3-velocity  is  v.  =  — —  , 

1  dt 

but  it  turns  out  that  in  relativity  theory  a  more  natural  choice 

for  the  relative  3-velocity  is  gotten  by  considering  the  4-vector 
dx 

A^  =  ,  and  then  considering  the  spatial  components,  A^  (i=l,2,3), 

of  A  .  It  can  be  shown  (see  Synge  -  "Relativity:  The  Special  Theory" 

y 

p.  130)  that 


v. 

1 

a 


a  +  ^)1/2 


Thus  letting  g  =  A  ,  v  =  v  we  see  that 

p  P 


v 

c 


a+V1/2 


To  relate  a(g,0)  to  the  conventional  cross  section  a(v, 0)  ,  we  choose 
the  time  axis  to  be  along  nd  ,  and  thus  in  the  center  of  mass  frame, 
ny  and  are  parallel.  Then  it  follows  from  1-19  that 


and 


dll  =  d  dP^  dP^  =  d  P  =  m  cosh  x  dw 


d'fi  =  d3’P  =  m  cosh  x  d'w 


But 


cosh  x  =  (1 


v 


4c 


.  (l+ig2)1/2 


Therefore  expression  2-41  becomes 

g  a (g,0)  N’N  (1  -  ~)1/2  d3P  (1  -  d3,P  dx 

4c  4c 

a (g , 0)  (1  +r  g2)~1/2  ’N  N  d3p  d3?p  dT  * 


v 
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Thus  we  see  that 


0( g,0)  (1  +  ~  g2)  1/2  =  a(v,9) 


so  that  we  get  the  conventional  expression  for  2-41,  i.e, 


v  a (v, 0)  NrN  d3PP3  P’  dx  . 


Chapter  III 

(1)  For  the  more  subtle  statistical  notions  see; 

D.  Ter  Haar,  "Foundations  of  Statistical  Mechanics",  Reviews  of 
Modern  Physics,  Vol.  27,  No.  3,  Section  A. 


Chapter  IV 

(1)  For  example  u^  could  be  defined  in  terms  of  the  particle  stream, 
i.e.  u^  parallel  to  M^  ,  or  in  terms  of  the  energy  flow. 

(2)  There  are  many  time  scales  involved  in  a  gas.  There  is  a  Kinetic 

time  scale,  x  ,  a  Hydrodynamic  time  scale,  x  .  The  Kinetic  time 
K  H 

scale  tells  us  how  long  the  gas  takes  to  get  to  a  local  equilibrium. 

The  Hydrodynamic  time  scale  gives  us  a  characteristic  time  for  a 

global  smoothing  out  by  Hydrodynamic  effects.  We  will  be  concerned 

with  times  in  the  interval  r  <  x  <  xu  ,  i.e.  before  x..  Hydrody- 

K  h  K. 


namics  is  not  applicable. 
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(3)  Eckart  -  Physical  Review,  58,  p.  919  (1940). 

(4)  Weinberg  -  Astrophysical  Journal  Vol.  168,  No.  2  page  1,  Sept.,  1971 

(5)  For  the  non  relativistic  discussion  of  this  see,  Landau  and  Lifshitz 
"Fluid  Mechanics"  (Addison-Wesley  1959). 

For  the  relativistic  discussion  of  this  see,  Weinberg  Astrophysical 
Journal,  Vol.  168,  No.  2,  Page  1,  Sept.  1971. 

(6)  See  Weinberg,  Ap.  J.  Vol.  168,  No.  2,  page  1,  Sept.  1971. 

(7)  See  Israel,  Journal  of  Mathematical  Physics,  Vol.  4,  No.  9,  p.  1163, 
Sept.  1963. 

To  see  that  the  two  definitions  are  equivalent  see  Appendix  B, 

(8)  See  Israel,  Journal  of  Mathematical  Physics,  Vol,  4,  No.  9,  Sept. 
1963,  equation  9-12. 

(9)  Ibid,  equation  9-13. 

For  the  explicit  calculation  see  Appendix  C. 

(10)  For  classical  discussion  of  Grad  Method  see  Grad,  Hb.  d.  Phys. 
Volume  12,  S.  Flugge  ed .  (Springer,  Berlin  1958). 

(11)  Anderson,  "Proceedings  of  the  Relativity  Conference  in  the  Midwest", 
Edited  by  Carmeli,  Fickler,  and  Witten  (Plenun  Press  1970),  page  109 


,  ■>  -i  r 
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(12)  Ibid,  see  page  117. 


(13)  Ibid,  see  page  119. 


(14)  See  Appendix  D. 


(15)  Israel,  J.  of  Math.  Phys. ,  Vol.  4,  No.  9,  Sept.  1963, 
equation  6.28. 


(16)  Ibid,  equations  6.35,  and  6.36. 

(17)  Ibid,  equation  8.14. 

(18)  Ibid,  equation  8.15. 

(19)  Ibid,  equations  6.29,  and  8.16. 

(20)  Consider  arbitrary  b°^  ,  let  b  =  ga^  b  then 

ag 


,ag  ,.ag  1  ag,  N  ,1  ag  , 
b  =  (b  -^-gb)+^-g  b 


BaB  +  |  gaB  b 


It  follows  then  that 


baB  P  P  +  C  =  Ba3(P  P  )  +  j  b(qa3P  P  )  +  C  = 
a  g  5  ag  4  ^  a  g  5 

=  Ba3P  P  +  C’  , 
a  g  5  * 

where  C*  =  \  b(-m^)  +  C-  and  g°^B  =  0 
5  4  5  ag 

Thus  we  can  absorb  the  trace  of  b  „  into  our  constant,  so  that 

ag 

we  can  without  loss  of  generality  let  b^g  =  0  • 


. 

•* 
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(21)  See  Appendix  E. 

(22)  Israel,  J.  of  Math.  Phys. ,  Vol.  4,  No.  9,  Sept.  1963,  equation  6.37. 

(23)  Ibid,  equations  6.35,  and  6.36, 

(24)  Ibid,  equation  6.37. 

Chapter  V 

(1)  For  detailed  calculations  see  Appendix  F. 

(2)  For  more  discussion  on  this  matter  see; 

Shucking,  Spiegel,  "Comments  on  Astrophysics  and  Space  Physics", 

Vol,  2,  p.  121  (1970). 

Stewart,  MacCallum,  and  Sciama,  "Comments  on  Astrophysics  and 
Space  Physics",  Vol.  2,  p.  206. 

Anderson,  "Proceedings  of  the  Relativity  Conference  in  the  Midwest" 
published  by  Plenum  Press,  1970,  page  109. 

(3)  Peebles,  R.A.S.C.  Journal,  Vol.  63,  No.  1,  page  5. 

(4)  Weinberg,  Ap.  Journal,  Vol.  168,  No.  2,  page  1,  Sept.  1971. 

(5)  Ibid,  see  equation  3.14. 
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(6)  Ibid,  see  equation  3.17. 

(7)  Ibid. 

(8)  Ibid. 

(9)  Ibid,  see  equation  3-19. 

(10)  Shucking,  and  Spiegel,  "Comments  on  Astrophysics  and  Space 
Physics",  Vol .  2,  p.  121  (1970). 
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APPENDIX  A 


We've  already  seen  that 

A-l  v  =  N  ds  dfl  , 

where  ds  and  dip  are  orthogonal  to  some  timelike  unit  vector  n 
If  we  choose  n  to  be  along  the  time  axis  then 

A- 2  ds  =  dx^  dx^  dx^ 

A- 3  dfi  =  d  Px  d  P2  d  P3  . 

Equation  A-3  is  a  volume  element  in  the  momentum  space, 
we  wish  to  deal  in  velocity  space  then  the  volume  element  is 


A-4  dll  =  d  u^  d  u^  d 

Then  it  can  be  shown  [see  J.L.  Synge  "The  Relativistic  Gas",  p.  16 
North  Holland  Pub.]  that 

3  5  2 

A-5  d£!  =  m  ■■■jjp  d U  ,  where  y  =  (1  -  ^y) 

c  c 


Eauation  A~1  now  becomes 

X 


A- 6 


N  -  y^  dx  dx  dx  du  du0  du. 


If 


v 
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Equation  A- 6  tells  the  number  of  particles  found  in  a  box  dx^  dx^  dx^ 
of  the  observers  space  at  time  t  ,  with  velocities  in  the  range 
du^  du^  du^  . 

5 

As  can  be  seen  that  because  of  the  y  term  the  distribution 
is  weighted  towards  the  high  speeds.  This  rather  "unnatural"  result 
disappears  if  we  stay  in  Momentum  space. 

Another  reason  why  velocity  space  is  unfavourable  is  because 
if  we  ever  want  to  treat  photons  it  couldn't  be  done  in  velocity  space, 
because  they  all  have  the  same  speed  c  and  therefore  there  are  no 
ranges  to  speak  of. 


Starting  with 
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APPENDIX  B 


B-l 


(f  =  x  Ayve 


,  v 


and  using  a  result  given  by  Israel  [J.  of  Math.  Phys.  Vol.  4,  No.  9, 
Sept.  1963,  equation  6,29] 


B-2 


a  =  M  o  ,  MS  1 

,v  2  ,v  4  f  P  v 
me  me  (pH - y) 


>  v 


me 


But  3  =  t~—  therefore  3 

k  T  ,  v 


me 
kT2  ,V 


T  so  that  equation  B-2  becomes 


B-3 


=  _  -0-  x  +  _JL_ 

,v  t2  ,v  2 

1  cl 


(y  +  -p 

C 


>V 


But  we've  already  seen  that  (see  equation  3-65) 


B-4 


u  =  - 

y 


(U  +  “~2> 


A  P 

p .  y  >v 


Substituting  equation  B-4  into  equation  B-3  and  then  substituting  equation 
B-3  into  equation  B-l  we  get 


QM  =  xamv  (- 


n  * 


rp  2  ,v  2  v 
T  c  T 


u.  ) 


which  simplifies  to 
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-  -  (% 

T 


(T 


T  u  ) 
2  v 


or 


B-5 


(T  +  T  u  ) 
,v  v 


Equation  B-5  is  just  the  phenomenological  form  of  the  heat  flux  vector 
where 
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APPENDIX  C 


We  know  that 


C-l 


oX 


V 


N  Y(P)  PV  dw  =  -  t 
o  4 


t  » 


W  N  fN 


I  6  (f)  5  (M7)  dco  d’co  dco 
o 


d  'co 


Let 


then 


\jr  _  p^-ph 


C-2 


3x 


v 


N  PXPUPVdw  =  -  y 
o  4 


r 


W  N  fN  <5(f)6(P  PU)  do)  d  *  go  dco  d'co 
o  o 


ot  8 

Substituting  f  =  P  PP  into  equation  C-2  we  get 


c-3  a 


Xyv 

o  ,v 


1  3 

-r  me 
4 


W  N  'N  <5  (fT,  )  6(P  Py  )  dio  d'u  dm  d  fco 
o  o  Horn 


7  b  « 
4  aS 


r  r  r 


W  N  'N  6(Pape)5(PXPU)  dco  d  *  a)  dco  d  * 


o  o 


*  .  * 
co 


But  6  ( f_,  )  =  0  ,  therefore 

Horn 


u  ,  =  -  b  .  Xa(Ul1  where 


O  V 


a$ 


aSAy  1  3 

X  =  -£■  me 


W  N  'N  6(Pape)£(PXPy)  dco  d'co  dco*  d'co* 


o  o 


. 
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APPENDIX  D 


For  the  sake  of  calculation  let  c  =  1  .  To  first  order  in 


D-l 


-M  Ma  =  -  (M  +  M  )  (M  a  +  M,  a) 
a  oa  la  o  1 


=  -  M  M  a  -  11%.  -  M-  M  a  =  p2 

oa  o  1  oa  la  o 


Qy  = 


i  (T0AV  +  T  Xv)(M  +  M  )[«>*  +  (MoX  +  Mu)(MoX  +  M  >*)] 

P 


-  -  (T  AV  +  T1XV)(M  +  M  )[<$y  +  \  (M  .M  U  +  M1XM  U  +  M  M  y)  ] 

p  o  1  ov  lv  A  2  oA  o  1A  o  oA  o 

P 

— — {(T  XVM  +  T  AVM,  +  T  AVM  )  [6y  +  ~  (M  XM  y  +  M  M  y  +  M  M1 

p  o  ov  o  lv  1  ov  A  2  oA  o  1A  o  oA  J 

P 

-  —  { (T  XvM  +  T  AvM  +  T  yvM  )  +  (T  AvM  M  My  + 
p  o  ov  o  lv  1  ov  2  o  oA  ov  o 

P 

+  T  AVM,  M  M  P  +  T  AVM  M  M  y  +  T  XVM  M  yM.,  + 

o  lv  oA  o  1  ov  oA  o  o  ov  o  1A 


+  T  XVM  M  M  y)} 
o  ov  oA  1 


yv  _  „  „y,  v 


*yv 


But  T  r =  y  uHuv  +  P  A'"  and  M^y  =  p  uy  ,  so  that  T^yVM^  %  =  -  p  yu 


y 


o 


o  ov 


T  yV  M  =  P  AyV  =  PM  y  ,  T  yVM  =  p  T  UVu  ,  M  yT  AVM  M  =  p3  y  uU 
o  lv  ly  1  1  ov  1  v  o  o  ov  oA 

T  AV  M  M  y  =  P  M,AM  M  y  =  0  ,  and  T  AvM  M  M,y  =  p2  y  M.U 

o  lv  oA  o  1  oA  o  o  ov  oA  x  l 


v)]j 
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Qy  =  -  -  {(-pyuy  +  PMxy  +  pT^u^)  +  —  (p3PuV1  +  P2P  M^)} 
P  V  P 


QP  =  -  (P±P)  M  y  -  T/V  u  . 
p  1  1  v 


Noting  that 


y+P 

P 


since 


u  u 

v  p 


pv 


=  0  then 


u 

v 


u 

V 
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APPENDIX  E 


E-l 


A  Z 

3  o 


ot3yv  4 

v  =  me  - 

o  33  93093  96 

a  3  y  v 


where 


Similarly 


z  (a, 3  )  =  4m  a 
o  y 


Kx(3) 


9z  9z„  a  o  k9(3) 

o  o  3  3  ,  2  0v 

-  -  =  4tt  a  - —  p 


33  33  33 

v  v 


6' 


£  tB-nKn(B)]  - 


-  B'n  Kn+1  (B) 


32z 


33  93 
y  v 


9  3  3^  ,  33y  _  yv 

33  3  ’  and  33  g 

y  v 


K„(3)  K9 (3) 

o  ,  3  nvny  ,  ,  2  yv 

=  4tt  a  - z—  3  3  +  4tt  a  - r—  g 


3' 


3 


E-2 


34z 


33  3 3 Q 9 3  93 
a  3  y  v 


K^)  rv  Q  li  K  (S) 

=  4tt  a  - =—  3  3  3  3^  +  4u  a  - — 


3' 


3 


[BVgaB  +  sVgUa  +  BVg'™  +  B°bV6  +  bV  VB 


K  ( g) 

„a„3  yv,  ,  ,  3  r  va  y3  ,  v3  ya  ,  a3  yv, 

+  3  3  g  ]  +  4tt  a  - r-  [gg  +gg  +ggJ 

3 
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Noting  the 

E-3 

E-4 

Equation  E 

E-5 


identities 


va  yg  ,  v g  y a  ag  yv  _  .va.yg  .vg.ya  ,  .ag.yv 
gg  +gg  +gg  =  A  A  +AAh+AAh 


yg  0  v  a  n  y  a  vg 

va  u  u  yg  u  u  vg  u  u  ya  u  u 

g  o  §  O  8  o  —  §  o 


v  y  . 
u  u  ag 


a  g 

u  u  yv 
— T~  § 


-  3 


a  g  y  v 

u  u  u  u 


v  y  ag  a  g  yv  ,  (v  y)  (a  g) 
uuA  +uu  +4uA  u 


6 


a  g  y  v 
u  u  u  u 


v  y  ag  ,  v  g  ya  ,  y  g  va  , 

=  u  u  g  +uug  +uug  + 


a  v  yg  ,  y  a  vg 

u  u  g  +  u  u  g 


,  a  g  yv 
+  u  u  g 


-2  simplifies  to 


Q  K,(g)  6K, (g)  3K  (g) 

agyv  ,  r  5  4  ,  3  n  a  g  y  v 

v  =  4tt  ma  { — - ~ —  H - ~ — >  u  u  u  u 

o  g  „3 


g 


g' 


2  fK4(3)  K3,  .  ( v  y) (a  g) 

+  4tt  mac  { - —  -  — 4  u  A  u  + 

g  g 

K  C 8) 

/  4 r  3V  r. va.yg  .  .vg.ya  ,  .ag.yv,  , 

+  4u  mac  { - x—  {A  A  +  A  A  +  A  A  }  + 


K  K 

,  /  2r  4  3-,  r  a  g.yv,  y  v.ag, 

f-  4tt  mac  {— r-  -  — {uuA  +  uuA  } 

g  g 


. 


. 
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APPENDIX  F 


We  know  that 


F-l 


,Pri^N2  1 

K  =  aF 


ft  =  3y  -  5  + 


3y 

nf3 


n  = 


K, 


K, 


Let  us  consider  3  ~  0  ,  i.e.  the  ultrarelative  limit  (T  ->  °°) 
3  ~  0  we  know  that  to  2nd  order 


K3(3) 


8 


1 

3 


8-3‘ 


F-2 


K, 


1 

2 


4-3" 

232 


so  that  we  have 


n  = 


K, 

YL, 


8-3  23i 


3' 


4-3' 


2  2  1 

f  (8-B  )  |(1 


-h  <8-s2)d  +  V 


3  2 


But 


_1_  =  _  32  dn 

Y-l  d3 


3 


4.1 
2  2  Y 


3  18 


for 


3  ~  0 


So  that  for  g 


0 


ft  =  3y  -  5  + 


3y_ 

nB 


,  but  nB  ~  4  + 


ft 


ft 


rv 

<"*w» 


2  3(—  +  -S— ) 

4  R  ^3  18; 

3(t  "  fs-)  ~  5  +  — - ^ 

4  +  3/2 
2 


k3  18' 

.2 


3' 


thus 


32\“1 


4  -  C  .  5  +  (4  +  C.)  1  (1  +  i^) 


F-4 


ft 


-1  - 


4 + j  (4 + 4> (i  -  4> 


ft  ~ 


2  02 

-1  +  — - 1-  1  - 

1  6  12 


ft 


JL_ 

12 


For 


we  have 
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It  should  be  noted  that  we  were  forced  to  take  to  2nd  order 

(i.e.  up  to  the  2nd  term  in  series  expansion)  in  calculating  Q  because 
of  the  factor  3  which  was  multiplying  p  in  this  expression  for  p  . 


We  also  know  that 


for  g  ~  0 


2  K2  K3  2  K3 

Pp  =  4 it  mac  —  —  =  4 it  mac  —  , 

3 2  K2  32 


F-5 


Pp  ~  32 7T  mac^  (-pr) 

0 


Let  us  next  calculate  A  for  3  ~  0  : 


F-6 


A  = 


3  2-1 

8tt  ma  c 


JJ 


3  ^3^y)  7 

sin  0  x —  g  a(g,0)  dg  d0 


where 


We  know  that 


y  =  26  (1  +  i  g2)172 


a(g,0)  =  a(v,6)  (1  +  -|  g2)1^2 


or 


,  .x  ya(v,0) 

o(g,e)  =  ^ — 


so  that  equation  F-5  becomes 


F-7 


A  = 


3  2  -1~ 

8tt  ma  c  a 

232 


3  7  ^3^^ 

sin  0  a(v,0)  d0  g  -  dg 


’ 


m 


-  136  - 


Near  3  ~  0  that  is  T  ->  °°  we  will  assume  that 


a(v,0)  =  a(0) 


so  that  equation  F-6  becomes 


F-8 


A  = 


o  3  2  -i 

ott  mg  c 

262 


K3(y)  ? 


g  d§ 


where 


a  = 


sin  0  a(0)  d0 


Consider  the  integral. 


F-9 


W  = 


K  (y) 

— -  g  dg 

•to  y 


us 


1  2  1/2 

ing  y  -  28  (1  +  ^7  g  )  we  transform  the  variable  of  integration 


to  y  i.e. 


2  tn2„  ,12,  2  y2  -  482 

y  =  48  (1  +  j  g  )  - >  g  =  1 - ~ - - 

8 

6  4  2  2  4  6, 

_ >  6  _  (y  -  12y  8  +  48y  8  ~  643  ) 

g 

66 

and  also  that  =  g  dg  ,  thus  equation  F-8  becomes 

6 


W  = 


+ 


88' 


48 


CO 

y6  K  (y)dy  - 

12 

6 

4 

y 

26 

6  j 

26 

roo 

y2K  (y)  dy 

_  64 

2  6 

3 

26 

K3(y)  dy  . 
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Consider  W  terra  by  term  we  have: 


,00 


I  = 


8 


fOO 


28 


28 


Y  K3(y) 


,8 


0 


K3(y)  y  dy 


y 


8 8  •'o 


/oo  b 

q  y  Kg  dy  is  bounded  so  that  we 


I  = 


c  2  8 

_1  _  J_ 

8  8 
8  8  j0 


K 


3  y  dy 


But  in  the  range  (0,28) 


thus 


K3(y)  , 


F-ll 


I  = 


~L  _  _8_ 

68'  68 


26 


0 


3  1  1 

y  dy  =  —  +  0(— ) 


8 


8 


8 


II  =  - 


12 

6  “3 

8  J28 


V  4  a  12 

KQy  dy  =  -  — 


12 


0 

■28 

0 


y4  k. 


k3  y 


By  similar  arguments  we  get 


F-12 


2  1 

II  =  -4  +  0  (  t)  , 


8 


8 


where  Cg  =  -12 


0 


y  Kg  dy  ,  For 


m-M 

B4  ;2B 


y  K3  dy  , 


K3  dy  , 

can  write 


For 

(y)  dy 

dy  • 


we  can  write 
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III  = 


48 


2  v  ,  48 

y  K  dy  =  — 

28  8 


,  2  8.  ,  48 

(y  K  -  — )  dy  +  — 

28  7  8 


28 


8  A 
-  dy 


which  can  also  be  rewritten  as 


III  = 


48 


0 


(  1  V  48 

(y  K  -  -)dy  -  — 

y  8 


28 


0 


(y2  K3  -  ^)  dy  + 


-  dy  . 
28y 


r°°  2  8 

The  integral  ~  J  q  ^  -  — )  dy  is  bounded,  and  in  the 

g 

interval  (0,28)  we  once  again  can  approximate  K^Cy)  ~  —  ,  therefore 

y 

integral  III  becomes 


F-13 


3  1 

III  =  -4  +  0(~r) 

6  8 


The  integral  IV  =  - 


64 


rco 


28 


dy  can  be  rewritten  as 


IV  =  - 


64 


,v  8  ,  U  A  64  /  8  1.  . 

(K3 - 2  +  dy - 2  (_T  ‘  ^  dy  > 

28  y  y  8  J28  y  y 


but  this  in  turn  can  be  rewritten  as 


IV  =  - 


64 

2 

8 


0 


(K  -  —  +  +  ~ 

K  3  2  yJ  2 

y  8 


28 


0 


(K3  ‘  "T  +  y5  dy 

y 


64 

62 


rco 


,8  ^  ^ 

(— 2 ' dy 
28  y  y 


By  arguments  similar  to  those  given  for  the  integral  II  we  get 


F-14 


IV  =  -  +  0(-^>  . 

8  8 


Using  the  results  given  in  equations  F-10,  11,  12,  13,  equation  F-9 


becomes 
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c.  c 

w  =  -£  +  0(-%  +  -|  +  oA-)  +  O(-i-)  +  0  A) 

B8  84  B6  B4  B4  B5 


Thus  we  finally  get 


w  =  Ur  +  0  A) 

B8  86 


or 


F-15 


W 


3 


Upon  substituting  equation  F-14  into  equation  F-7,  and  together  with 
equations  F-4,  5,  equation  F-l  becomes 


^  .16c  m,  n3 
k  »  ( - )  3 

9TTOC, 


F-16 


